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. ABSTRACT 

' We use Spitzer Space Telescope and Herschel Space Observatory far-infrared data along with 

\ ground-based optical and near-infrared data to understand how dust heating in the nearby 

. face-on spiral galaxies M81, M83, and NGC 2403 is affected by the starlight from all stars 

' and by the radiation from star forming regions. We find that 70/160 /im surface brightness 

ratios tend to be more strongly influenced by star forming regions. However, the 250/350 /im 
and 350/500 /im surface brightness ratios are more strongly affected by the light from the total 
stellar populations, suggesting that the dust emission at >250 /xm originates predominantly 
from a component that is colder than the dust seen at < 160 /im and that is relatively unaffected 
' by star formation activity. We conclude by discussing the implications of this for modelling 

H \ the spectral energy distributions of both nearby and more distant galaxies and for using far- 

^ infrared dust emission to trace star formation. 

Key words: galaxies: ISM, galaxies: spiral, infrared: galaxies, galaxies: individual: M81, 
galaxies: individual: M83, galaxies: individual: NGC 2403, 



1 INTRODUCTION 

Since the launch of the Infrared Astronomical Satellite (IRAS; 
iNeugebauer et al.l Il984h . astronomers have frequently used 
far-infrared emission to trace star formation in extragalactic 
sources. Far-infrared emission has been very popular to use 
for estimating star formation rates because it is not as strongly 
affected by dust extinction as ultraviolet and optical tracers of 
star formation and because many surveys, particularly surveys 
of high-redshift sources, have detected a multitude of galaxies in 



the far-infrared. Several equations have been derived to calculate 
star formation rates either from far-infrared measurements by 
themselves or from far-infrared data com b ined with optical 
or ultraviolet da t a Jscoville & Younl Il983l : iBuat & Xul 1 19961: 
iKennicuttI Il998j: iKennicutt et al.l l2009h . While it is commonly 
assumed in the use of these equations that dust is primarily if not 
solely heated by star forming regions, the heating source for the 
dust producing far-infrared emission has been widely debated for 
decades. Following the completion of the IRAS surveys, several 
papers showed that IRAS 60 and 100 /im emission was strongly 
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corre la ted with star fomatio n activity (e . g . Devereux & Youn^ 
119901 ; iDevereuxetai] 1 19951 : iBuat & Xul 119961) , while others 
argued that a significant fraction of dust emission originated 
from d ust heated by evolved stars (e.glLonsdale Persson & Heloul 
19871; IWalterbos & Schwerind 1 19871: ISauvage & Thuaij 1 19921 : 



Walterbos & Greenawaltl Il996l : iKong et al.l |2004) . Add ition 



ally, analyses per formed by iDevereux & Youn J ( Il992h and 
IPevereux & Youn^ (1993) on IRAS data combined with 160 or 
170 /xm data from the Kuiper Airborne Observatory indicated 
that emission at >100 fxm also originated from dust heated by 
star formation. Later observati ons of nearby spira l galaxies with 
the Spitzer Space Telescope jWemer et alj |2004|) demonstrated 
that optical and ultraviolet star formation tracers were very 
strongly correlated with 24 ^tm dust emission I calzetti et al. 2005, 



2007 : iPrescott et al.l 120071 : iKennicutt et al. 12007^ .2009,: .Zhu et al.l 
2008 ) or with total infrared emission Kennicutt et al.l l2009h . 



However, similar analyses showed that the correlation between 
optical/ultraviolet star f ormation tracers and either 70 or 160 pim 
emission was weaker jCalzetti et aP 1201 Oh . which implied that 
not all of the far-infrared emiss ion origina t ed fro m dust heated 
by star formation. Additionally, iHinz et 2A\ ( l2004h demonstrated 
that large scale 160 /^m emission from M33 qualitatively appeared 
more similar to large scale K-band emission than to large scale 
Hq, 24 pm, or 70 /xm emission, implying that the dust emitting at 
160 was at least partly heated by the evolved stellar population. 
Furthermore, some dust models ap plied to IRAS, Infrared Space 
Observatory ( 'Kessler et al.l Il996h , Spitzer, and ground-based 
submillimetre data for nearby galaxies have suggested that heating 
by evolved stars was needed to reproduce th e observed spectral 
energy distributions for these galaxies (e.g. Draine et alj |2007l : 
IdaCu nha et al. 2008; Bianchi 2008; Popsecu et al. 201 1). 

The first resu lts from the Herschel Space Observatory 
jPilbratt et al.l201(]h brought attention back to this issue mainly be- 
cause the telescope provided data at longer wavelengths than what 
had been provided by either IRAS or Spitzer and at a higher signal- 
to-n oise than what was p ossible with ground-based instrumenta- 
tion. iBendo et alj feOlObf) demonstrated that colours between 160 
and 500 fim in M81, an Sab galaxy ide Vaucouleurs et ^Il99lh . 
were more strongly dependent on galactocentric radius than on far- 
infrared surface brightness, which implies that the dust was heated 
primarily by the evolved stars in the galaxy (which have a sur- 
face brightness that varies primarily with radius) rather than by star 
forming regions (wh ich are primari ly found in the infrared-brigh t 
regions). In contrast, iBoquien et alj t201Q) and lVerlet et al.l ( [20101) 
found that 1 00-250 /xm emission from c ompact sources in M33, an 
Scd galaxy dde Vaucouleurs et al. |[l99l|), were strongly correlated 
with star formation as measured using Ha emission, 24 fim emis- 
sion, or a combination of the two. 

The goal of this paper is to follow-up these early Herschel 
results using three nearby face-on spiral galaxies with differing 
Hubble types taken from the Very Nearby Galaxies Survey, a Her- 
schel guaranteed-time photometric and spectroscopic survey of 13 
galaxie s with t he Photodetector Array Camera and Spectrometer 
(PACS: |Pogiitsch et al 2010) and Spectr al and Photometric Imag- 
ing Receiver (SPIRE: lOriffin et aLllioiOh . We will primarily com- 
pare surface brightness ratios between two Herschel bands to both 
stellar surface brightness and to star formation rate. Surface bright- 
nesses measured in single bands may be correlated to star formation 
through dust heating, in which the dust is heated by the star forming 
regions and becomes more lumin ous, or through the Schmidt law 
( ISchmidlll95'^ : lKennicutlll998b '), in which star formation is corre- 
lated to the amount of dust present. In contrast, surface brightness 



ratios are advantageous to use in that they are largely independent 
of the total dust surface density. Instead, surface brightness ratios 
will primarily depend only on the average temperature of the dust 
emitting at those wavelengths, with temperature variations caused 
either by variations in the amount of a warmer dust component rel- 
ative to a colder dust component or variations in the overall tem- 
perature of a single dust component. This is similar to the approach 
used by Bendo et al. (2010b), but the advantage of the analysis in 
this paper is that we will not rely upon galactocentric radius as a 
proxy for the surface brightness of the evolved stellar population 
and that we will use a better tracer of star formation. 

Section |2] describes the three galaxies in our sample (M81, 
M83, and NGC 2403) and why they are well-suited for this analy- 
sis. Section[3]provides an overview of the data from Herschel and 
from other sources as well as details on the preparation of the data 
for the analysis. The analysis itself is presented in Section|4l the im- 
plications of the results is discussed in Section|5] and a summary is 
provided in Section|6] 



2 SAMPLE GALAXIES 

As stated above, the three galaxies used for this analysis were all 
galaxies selected from the Very Nearby Galaxies Survey (PI: C. 
Wilson), a guaranteed-time survey of 13 well-studied nearby galax- 
ies with a diverse range of properties. The galaxies for the analysis 
in this paper were selected because they are all spiral galaxies with 
inclinations of <60° from face-on, they have optical discs larger 
than 10 arcmin, and they are at distances of less than 5 Mpc, so 
sub-kpc structures can be easily studied within the galaxies. The 
galaxies are also roughly representative of early-, mid-, and late 
type spiral galaxies. Details on the three galaxies used in this anal- 
ysis are given in Table[T] The images used in the analysis are shown 

in Figures [TTSl 

^MSl (NGC 3031) is an Sab galaxy ids Vaucouleurs et all 

Il99ll) that is th e brig htest galaxy in the M81 Group 
dKarachentsev et ^ 1 20021) . As this is an early-type spiral 
galaxy, the bulge is a relatively large and bright component of 
the galaxy. The dust emission is unusually extended compared to 
other early- type spiral galaxies; the dust dis tribution is similar to 
late-type spiral galaxies (Bendo et alj l2007h . The dust emission 
also does not qualitatively appear similar to the stellar emission. 
Furthermore, the star forming regions lie primarily in the outer 
disc, whereas the centre of the galaxy is relativel y devoid of star 
formation, as shown b y the Ha image (also see IPevereux et al.l 
1 19951 : 1 Allen et al]|l997i' ), although incompletely-subtracted contin- 
uum emission from the bulge is visible in the image. Consequently, 
this galaxy is very useful for differentiating between effects related 
to heating by starlight and effects related to heating by star forming 
regions. The gal axy does contain a low luminosity active g alactic 
nucleus (AGN; iMarkoff et ai]|2008l : iMoustakas et alJliolO ) that 
could affect the observed spectral energy distribution (SED) for 
the nuclear region. We comment on this in the analysis where 
appropriate. Also, the 250-500 /im images show that the galaxy 
is surrounded by extended e mission that was d etermined to be 
foreground cirrus emission by ('Pavie s et"al .120101) . The foreground 
structures are fainter than much of the emission from M81 itself, 
so by choosing appropriate surface brightness thresholds for our 
analysis, we avoid working with data that is strongly affected by 

the foreground emission. 

M83 (NGC 5236) is an Sc galaxy ('de Vaucouleurs 
with a grand design appearance (Elm egreen & Elmegr 



s etal."l99l|) 

green mt}. 
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Table 1. Basic data on sample galaxies. 



Galaxy 


Hubble 


Size of Optical 


Position Angle 


Inclination 


Distance 


Physical Size 




Type " 


Disc (arcmin)'' 


of Major Axis*"^ 




(Mpc) 


of 36 arcsec (pc)'* 


M81 (NGC 3031) 


SA(s)ab 


26.9 X 14.1 


157° 


59.0° 


3.6 ±0.4-'' 


630 ± 70 


M83 (NGC 5236) 


SAB(s)c 


12.1 X 11.5 


45° 9 


24° f 


4.5 ±0.2'* 


790 ± 30 


NGC 2403 


SAB(s)cd 


21.9 X 12.3 


127° 


62.9° <= 


3.2 ±0.3-'' 


560 ± 50 



" Hubble types are from lde Vaucouleurs et al. ( 1991). 

' Data are for the D25 isophote as given by de Vaucouleurs et alj jl99lb except where specified. 
Position angle is defined as degrees from north through east. 

This angular size represents the size of the bins used in the analysis in this paper. See Section [X4l for information on the 
selection of this bin size. 



^ Data are taken from'de BoketaD ^iooi). 
f Data are taken from Freedman et al. 1 200 lb . 

^ Data are taken from Comte 1 1981^^ 

Data are taken from iThim et al., >2003l) . 




Figure 1. The 30 x 20 arcmin images of M81 used in this analysis. The Hq emission traces the gas that is photoionised by star forming regions (although, in 
the case of M81, some residual continuum emission from the bulge is present). The 1.6 fim image traces the total stellar population. The 70-500 fim emission 
trace the thermal emission from colder dust within the galaxy. North is up and east is left in each panel. The green circles in the lower left corner of each panel 
shows the FWHM of the data. The cyan ellipse in the 1.6 /im image shows the optical disk of the galaxy. 
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Figure 2. The 20 x 20 arcmin images of M83 used in tliis analysis. See tlie caption of Figure[T]for other information. 



While the centre is the site of a strong starburst (e.g . iTalbot et al 



19791 : iTelesc o & Han^eil Il98d: ICondon et alj Il982l : iBohlin et al.. 
1983 : Trinchi eri et al.1 Il985l: iTumer et al.1 Il987h . star formation 



rates in the disc are not as extreme. Unlike the other two galax- 
ies in this analysis, M83 is a case where the stellar emission and 
dust emission look qualitatively similar. However, it should be rep- 
resentative of many other nearby grand-design Sb-Sc galaxies. 

NGC 2403 is an Scd galaxy Jde V aucoul eurs et alj|l99lh wit h 
a flocculent appearance telmegreen & Elmegreeij 1 19821. Il987h . 
The galaxy is the secon d brightest galaxy in the M81 Group 
jKarachentsev et aj]|2002[) . The extended, asymmetri c dust emis- 
sion in this galaxy is typical of many Sc-Sd galaxies iSendo et al.l 
[2007.) . The distribution of star forming regions in NGC 2403 differs 



significantly from the distribution of starlight; while the starlight 
peaks in the centre of the galaxy, the regions with the strongest star 
forming activity are actually seen off-centre, as seen in Ha shown 
here as well as Ha and 24 fim images shown in previous publica- 
tions (e.g. Drissen et al. 1999; Bendo et al. 2008). The star forming 
region with the highest Ha and 24 emission, which is l abeled 
as region 44 (VS 44) in the catalog ue oflVer on & Sa uvayrd il965l) 
and region 128 in the catalogue bv lHodg~ & Kennicujt 1 1983h . is 
located northeast of the centre of the galaxy. Since this galaxy has 
multiple bright, off-centre star forming regions, it is extremely use- 
ful for disentangling effects dependent on star formation from ef- 
fects dependent on either radius or starlight. 
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Figure 3. The 21 x 18 arcmin images of NGC 2403 used in this analyis. The stai' forming region with the highest Ha intensity (VS 44) is marked with a green 
square. See the caption of Figure [T]f or additional information. 



3 OBSERVATIONS AND DATA REDUCTION 

3.1 Far-infrared and submillimetre data 

3.1.1 MIPS 70 ^im data 

Although 70 /im images were produced at the same time as 160 /im 
images with PACS, we decided to use MIPS 70 data for M81 
and NGC 2403. PACS data have a smaller point spread function 
(PSF) than the MIPS data, so PACS data can be used to map struc- 
tures on finer scales. As we will be matching the PSFs of the 70 
images to the wider PSF of the SPIRE 500 pim images, the PSF size 
is not important for our analysis. Comparisons between the PACS 
and MIPS 70 fim data for our sample showed that the MIPS data 
are more sensitive to extended emission, which is more important 
for our analysis (although both the PACS and MIPS 160 fim im- 
ages have similar sensitivities). However, the MIPS 70 data 
suffer from latent image effects from bright sources that manifest 
themselves as dark streaks in the MIPS data. Hence, the PACS data 



are preferable to use for brighter sources, while the MIPS data are 
preferable to use for fainter sources. In the M81 and NGC 2403 
data, the latent image effects are less problematic than the lower 
sensitivity of the PACS data, which is why we decided to use the 
MIPS data for the analysis on those galaxies. In images of M83, 
however, the PACS data have a sufficiently high signal-to-noise 
that the sensitivity issues are not as much of a concern, and the 
MIPS images are strongly affected by latent image effects, so we 
use PACS 70 fim data for that galaxy. 

The 70 fim images for M8 1 and NGC 2403 were a cquired with 
the M ultiband Imaging Photometer for Spitzer (MIPS; iRieke et al] 
l2004h a s part of the Spitzer Infrared Nearby Galaxies Survey 
(SINGS: lKennicutt et al.ll2003h . These observations consist of two 
scan map observations performed using the medium scan rate (6.5 
arc sec s^^). The raw data from the Spitzer archive were repro- 
cessed using the MIPS Data Analysis Tools dCordon et"aLll2005h 
along with additional processing steps. First, ramps were fit to the 
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reads to derive slopes. In this step, readout jumps and cosmic ray 
hits were also removed, and an electronic nonlinearity correction 
was applied. Next, the stim flash frames (frames of data in which a 
calibration light source was flashed at the detectors) were used as 
responsivity corrections. After this, the dark current was subtracted 
from the data, and an illumination correction was applied. Follow- 
ing this, short term variations in the signal (often referred to as drift) 
were removed, and additional periodic variations in the background 
related to the stim flash cycle were subtracted; this also subtracted 
the background from the data. Next, a robust statistical analysis was 
applied to cospatial pixels from different frames in which statisti- 
cal outliers (which could be pixels affected by cosmic rays) were 
masked out. Once this was done, final mosaics were made using 
pixel sizes of 4.5 arcsec. The residual backgrounds in the data were 
measured in regions outside the optical discs of the galaxies and 
subtracted, and then flux calibr ation factors (given as 702 MJy sr~^ 
[MIPS instrumental unit]^^ by [Gordon et alj fzOOTi) ) were applied 
to the data. An additional nonlinearity correction given as 

/70Mm(true) = 0.581(/70M™(measured))^'^^ (1) 

bv lDaleetal] ( l2005l) was applied where the surface brightness ex- 
ceeded 66 MIy sr~^. We also applied a colour correc tion of 0.901 
for a 30 K blackbody given by IStansberrv et akl ( 12007,) : the actual 
dust temperatures may differ from this value, but we anticipate that 
the associated colour correction will be within 10% of this value. 
The final data have a point spread function (PSF) wit h a full-width 
half-m aximum (FWHM) of 18 arcsec according to [Gordon et al.l 
( |2007[) . Calibration and non linearity uncertainties are estimated to 
be 10% l lGordon et al]|2007l) . 

3.1.2 PACS 70-160 fim data 

The PACS observations were performed as pairs of orthogonal 
scans using a 20 arcsec s~^ scan rate. PACS can perform simul- 
taneous observations in only two wave bands; the VNGS chose the 
70 and 160 fim bands since they were expected to bracket the peak 
of the SED better. Four pairs of observations were performed on 
each galaxy. The observations covered areas of 40 x 40 arcmin 
around M81 and NGC 2403 and 25 x 25 arcmin around M83 so 
as to sample an area with a size that is at least 1.5 times the optical 
disc. 

The PACS data were reduced with the H erschel Int eractive 
Processing Environment (HIPE) version 5.0 l2O10h and a 
pipeline adapted from the official one. After applying the standard 
tasks that reformat the telemetry, associate them with the correct 
pointing information, and flat-field and calibrate the signal into Jy 
pixeP^, we have applied the following specific treatments. First we 
deal with the possibility for electrical crosstalk, which potentially 
affects column of all matrices (groups of 16 x 16 bolometers 
in the PACS arrays) when a bright source illuminates column 15, 
by systematically masking columns 0. Then we proceed to scan 
the signal for glitches (cosmic rays) and other possible outliers. 
We use the so-called second-level method, where what we scan 
are not the individual pixel timelines but rather reconstructed time- 
lines made of all samples acquired on a given position of the sky. 
We use the timeordered option that is more robust to the pres- 
ence of strong brightness gradients in the sky. We check that we do 
not over-deglitch the data and do not flag by mistake object struc- 
tures as glitches by projecting the glitch mask in the sky. In this 
glitch "map" the distribution of glitches should be discontinuous 
and the object structure should not be recognisable. A satisfactory 
setting for the na parameter of the task is 25-30 which leads to 



a glitch rate of less than 1% of the data. We produce final maps 
using Scanamorpho43 (Roussel, submitted) with its standard 
settings, as Scanamorphos was the best mapmaking software avail- 
able to us for removing low-frequency drift and reproducing ex- 
tended emission in PACS data. After this, the median residual back- 
grounds were measured in ~ 1-4 arcmin wide strips near the edges 
of the observed regions and outside of the optical discs. These back- 
grounds were then subtracted from the data. 

In each map, north is up, east is left. The pixel size was set 
to the default value of 1.4 arcsec at 70 fim and 2.85 arcsec at 
160 /xm, which corresponds to 0.25 times the FWHM of the PSF 
in typical PACS observations. The PSF has a 3-lobed structure and 
pronounced Airy rings; the FWHM of the major axis in individ- 
ual scan maps in these data is 6 arcsec at 70 /im and 12 arc- 
sec at_j_60_jini_according_tojhe PACS Observer's Manual from 
the lHerschel Space Observatorvl l l2010a|{'^ The calibration uncer - 
tainty is 10% at 70 ^m and 20% at 160 ^Im jPogUtsch eral]|2010l) . 
Colour corrections and other effects are expected to be small com- 
pared to the calibration uncertainty and are thus not included here. 

3. 1.3 SPIRE 250-500 fim data 

The SPIRE observations for each galaxy were performed as a pair 
of orthogonal scans using a 30 arcsec s~^ scan rate and nominal 
bias voltage settings. The areas covered in the scan maps, again 
selected to sample an area at least 1.5 times the size of the optical 
disc, are 40 x 40 arcmin for M81, 20 x 20 arcmin for M83, 30 x 
30 arcmin for NGC 2403. 

The SPIRE 250-500 /im observations produce timeline data 
that were processed using a customized versi on o f the official 
scan m ap pipeline script (see lGriffin et al.l ( l2009h and lDowell et al.l 
( I2OIO ) for more information) run in a version of HIPE with the 
continuous integration build number 4.0.1343, which was the de- 
velopers' branch of the data reduction software. This version of the 
pipeline is similar to the script in the publicly-released version 5 of 
HIPE, and we also used the new flux calibration product (based on 
Neptune observations) that was included in version 5. 

The first steps of the timeline processing apply concurrent 
glitch removal (the removal of cosmic rays that affect all detec- 
tors in an individual array) and then wavelet glitch removal (the 
removal of cosmic rays from individual detectors), with the mod- 
ule adjusted to mask 7 samples following a glitch. Next, we applied 
an electrical low pass filter correction and flux calibration, and then 
we re-applied the wavelet glitch removal, as, in the version of the 
software that we were using, an additional deglitching step found 
to improve the removal of glitches from the data. After this, we 
applied a time response correction. 

We skipped applying the default temperature drift correction 
and baseline subtraction and instead used a custom method called 
BriGAdE (Smith et al., in preparation) to remove the temperature 
drift and bring all bolometers to the same level. For each bolometer 
in each array, BriGAdE fits a linear function to the relation between 
the bolometer signal timeline and the array's thermistor signal time- 
line after 'jumps' (sudden DC offset in the thermistor timelines) 
have been removed. The timeline fitting is applied to the bolome- 
ter data for the whole observation including turnaround data but 
excluding samples with signals from bright sources or unidentified 

Available at|http://www2.iap.fT/users /roussel/herschel/|. 

The PACS Observer's Manual is available at |http://herschel.esac.esa.iiit/ 1 
Docs/PACS/pdf/pacs_om.pdf . 
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glitches. For the 250 and 500 /im arrays, where we have two func- 
tional thermistors, we used the thermistor providing the best fit to 
each bolometer. The 350 /im array has only one functional ther- 
mistor, and so that thermistor needed to be used for this processing 
step. The resulting functions, which directly relate the thermistor 
signal to the background signal (including drift) measured by the 
bolometers, are then used to calculate and subtract the background 
signal from each bolometer signal timeline. We have found that this 
method improves the baseline subtraction significantly, especially 
in cases where there are strong temperature variations during the 
observation. 

We used the naive mapper in HIPE to create the final maps 
used in this analysis. We set the pixel size to those listed in Table|2] 
because they are the smallest pixel sizes possible that do not in- 
clude significant numbers of map pixels that were not crossed by 
detectors. We then subtracted the median background levels mea- 
sured in ~ 4-6 arcmin wide regions outside the optical discs of 
the galaxies near the edges of the mapped area. The SPIRE PSFs 
have FWHM that vary as a function o f pixel scale according to the 
SPIR E Observer's Manual from the iHerschel Space Observatory! 
feOlOb'PI: Table|2]lists the FWHM for the pixel sizes that we used. 
By default, the pipeline produces monochromatic flux densities for 
point sources where i/f^, is constant, so we needed to apply cor- 
rections that adjust the monochromatic flux density values from 
point source to extended source val ues and colour corrections (see 
IHerschel Space Observatory! 201 Obi , for more information), both of 
which are also listed in Table |2l The 350 fim data are also multi- 
plied by 1.0066 to correct for an issue related to an update of the 
filter profile, and the other corrections are based on the updated 
filter profile. The colour corrections that we use assume that the 
sources are extended and have dust SEDs that resemble blackbod- 
ies modified wi th emissivity funct ions that scale as (which is 
consistent with iLi & Draindl200lh . The temperatures are assumed 
to lie between 15 and 25 K. We compared the slopes of these func- 
tion at the central wavelengths for each SPIRE band to the slop es in 
the colo ur correction tables given by Herschel Space Observatory! 
( l2010bl) to determine the colour corrections given in Table |2| The 
final calibration uncertainties in the data include 5% systematic un- 
certainties that is correlated across all SPIRE bands and 2% rand om 
uncertainties in each band faerschel Space Observatoryl!2010bl) . 



3.2 Tracer of total stellar emission 

As a tracer of the surface brightness of the total stellar popula- 
tion, we used 1.6 ^lm (H-band) images from the "2MASS Large 
Galaxy Atlas l ljarrett et"al ]!2003l) . These data sample the Rayleigh- 
Jeans side of the thermal emission from the total stellar populations 
within galaxies, and they are relatively unaffected by dust extinc- 
tion. Shorter wavelength data would tend to be strongly affected by 
dust extinction effects and could be affected more by variations in 
the ages of the stellar populations. Longer wavelength data, particu- 
larly 3.6 fim data from Spitzer, is even less affected by dust extinc- 
tion than the 1.6 fim data. However, multiple recent studies have 
demonstrated that the 3.6 ^m band may sometimes include ther- 
mal emission from hot dus t at 700-1000 K associated with sites of 
very strong sta r formation ( Smith & Hancocdl2009l : !Mentuch et al.l 
l2009i,.2010.) . .Mentuch et al. (1010) even showed that such emission 
has been detected in the 3.6 /im images of M81 and NGC 2403. 



20 The SPIRE Observer's Manual is available at http://' 
herschel. esac.esa.int/Docs/SPIRE/pdf/spire_om.pdf . 



Furthermore, iFlagev et al.l ( |2006|) have shown that polycyclic aro- 
matic hydrocarbon spectral feature emission may also contribute 
up to 50% of the observed emission in the 3.6 nm Spitzer band. 
Hence, the 1.6 ^m band appears to be the best band to use for trac- 
ing total stellar populations, as it is the best compromise between 
dust extinction effects and dust emission effects. 

The processed, calibrated images were acquired from the 
NASA/IPAC Extragalactic Database. The 1.6 /im images all have 
pixel sizes of 1 arcsec pixel" Th e FWHM of the PSF is reported 
to be 2-3 arcsec djarrett et al.l2003l) : we use 2.5 arcsec for our anal - 
ysis. The flux calibration uncertainties are 3% djarrett et al.l!2003!) . 
To prepare the images for analysis, we measured median back- 
ground surface brightnesses in multiple circular regions outside the 
optical discs of each galaxy and subtracted from each image. We 
then identified bright foreground stars as unresolved sources with 
/i.6fim//24fim >50 (uslug rcprocessed 24 fim data from obser- 
vations originally performed by SINGS and by 'Eng elbracht et al] 
(l2005|)) and removed them by interpolating over them. 

3.3 Tracer of star formation 

We use Ha images of the galaxies as a straightforward tracer of 
how star formation heats the ISM in these galaxies. These images 
trace both the photoionised gas in star forming regions themselves 
as well as diffuse emission from gas that is ionised by photons es- 
caping from the star forming regions while tending not to include 
emission from other sources. However, the Ha images will be af- 
fected by extinction, with the extinction effects expected to be the 
strongest in star forming regions. This may flatten some of the 
relations between Ha intensity and other quantities. It would be 
preferable to use a tracer of star formation that corrects for dust 
extinction by combining the Ha emission (or emission in another 
band) with dust emission associated with star formation, such as 
24 fim emission. However, these techniques have only been de- 
veloped for either compact sources within galaxi es dCalzetti et al.l 
2007.) or for global flux density measurements dZhu et a" 20081 : 
Kennicutt et al.l2009l) . No corrections have been created for diffuse 
emission within galaxies, and using dust emission in any band to 
attempt this may be difficult, as even diffuse mid-infrared emission 
can poten tially contain diffuse dust emission h eated by evolved 
stars (e.g. lArendt et al.lll998l : lLi & Drainell200ll) as well as stellar 
emission. However, analyses with Spitzer data have shown that very 
few Ha sources in nearby spiral galaxies are severely obscured by 
dust (Prescott et al. 2007), and so even without dust extinction cor- 
rections, the Ha data should still perform reasonably well at tracing 
emission from star forming regions. Although we may expect the 
data to exhibit some additional scatter without extinction correc- 
tions, the uncorrected Ha images are still the most suitable data for 
tracing how star formation regions heat both compact and diffuse 
regions. 

The Ha images f o r M81 and NGC 2403 were originally made 
by iBoselU & Gavazzil d2002l) using a 1024 x 1024 CCD with a 
0.69 arcsec pixel" ^ scale at the 1.20 m Newton Telescope at the 
Observatoire de Haute Provence. The observations for each galaxy 
consist of a set of separate pointings that completely cover the op- 
tical disk of the galaxy. The observations were performed with in- 
terferometric filters where one was an in-band filter and the other 
was an off-band filter; the Ha image was made by subtracting one 
image from the other. The FWHM of the PSF was reported to vary 
between 2 and 4 arcse c; we use 3 arcsec. Calibr ation uncertainties 
in the data are 5%. See iBoselU & Gavazzil d2002l) for additional de- 
tails. 
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Table 2. Characteristics of the SPIRE 250-500 /jm data and corrections applied 
to the data. 



Wave Band 


Pixel 


PSF 


Correction from Point 


Colour 


(/im) 


Scale 


FWHM 


Source to Extended 


Correction'' 




(arcsec) 


(arcsec)" 


Source Emission'' 




250 


6 


18.2 


0.9939/1.0113 


0.993 ± 0.009 


350 


8 


24.5 


0.9920/1.0087 


1.000 ±0.007 


500 


12 


36.0 


0.9773/1.0065 


1.000 ± 0.008 



" The FWHM in the map data varies with pixel scale. These values are for the 
pixel scales that we used. 

^ These are multiplicative correc tions based on the values given by the 
iHerschel Space Observatory! )2010bl) . 



The Ha image for M83 was originally made bv lMeurer et alj 

1 I2OO6I) as part of the Sui^ey for Ionization in Neutral Gas Galaxies 
(SINGG) and was distributed as part of SINGG release 1.0. The 
observations were performed at the Cerro Tololo 1.5 Meter Tele- 
scope using the 2048 x 2048 CFCCD, which has a plate scale of 
0.43 arcsec pixel" ^. On- and off-band images were obtained using 
narrowband filters; the Ha image was made by subtracting the on- 
band image from the off-band image. The reported median FWHM 
of the PSF of the survey observations is 1.6 arcsec. C alibration un- 
certai nties are 4%. Additional details are provided bv lMeurer et al.l 
( |2006() . 

To prepare the Ha images for use in our analysis as a tracer of 
both the compact and diffuse emission from photoionised gas, we 
identified foreground stars in the image, which appeared as incom- 
pletely subtracted point sources, and also measured and subtracted 
the median background emission in circular regions at the edges 
of the images (in the case of M81 and M83) or at the edge of the 
region with uniform coverage (in the case of NGC 2404). Addition- 
ally, for M81, most of the emission within the central 5.7 x 3.0 ar- 
cmin (approximately the central 3 kpc) of the Ha image originates 
from either the AGN or incompletely subtracted stellar emission 
from the bulge; it even has a light profile with artefacts that are 
similar to the incompletely-subtracted foreground stars. This area 
is masked out for the analysis. We corrected the Ha measurements 
for foreground dust extinction within the Milky Way using calcula- 
tions from th e NASA/IPAC Extragalactic Database that are based 
on data from ISchlegel et al"] l[r998); the R-band extinction values 
are given in Table [3] We also applied corrections to remove the 
[Nil] line emission; we assume that the [Nll]/Ha ratios are approx- 
imately uniform across th e discs of the galaxies. The Ha data from 
iBoselh & Gavazzil ( l2002l) contain emission from both the 6548 and 
6583 A lines, while the data from Meurer et al. (2006) contain 
emission from only the 6 583 A line. In the case o f M81, we used 
data on HII regions from iGamett & Shieldsl ( 1 19871) to estimate the 
[Nil] (6548 and 6583 A)/Ha ratio to be 0.40 ± 0.13. For M83, the 
[Nil] (65 83 A)/Ha ratio is estimated to be 0.34 based on the data 
from Boissier et al] |2(K )5). From the radial strip data presented by 
iMoustakas et alj fcOlOl) . we calculated the [Nil] (6548 and 6583 
A)/Ha ratio for NGC 2403 to be 0.28 ± 0.05. 



3.4 Convolution and binning 

To compare images from different wave bands without worrying 
about beam effects, it is necessary to match the PSFs of all of the 
data to the PSF of the wave band with the lowest resolution. Oth- 
erwise, the analysis will be strongly affected by variations related 



Table 3. Foreground extinction values used to correct Ha fluxes. 



Galaxy 




M81 


0.214 


M83 


0.176 


NGC 2403 


0.077 



" These values were calculated by the NASA/IPAC Extragalactic Database 
using data from lSchlegel et alj h998l) . These extinction values are only for 
the foreground and do not include extinction within the target galaxies. 



to PSF shapes, with side lobes appearing around bright sources in 
maps based on surface brightness ratios and with fainter regions 
appearing biased towards emission from the wave bands with the 
broader PSFs. For our analysis, we matched the PSFs of the data to 
the PSF of the 500 ^m data, which has a FWHM of 36.0 arcsec in 
the maps. 

The SPIRE PSFs are very close to Gaussian (i.e. the Airy rings 
are weak compared to the peak of the PSF), so we could use the 
equation 



K{r) 



(2) 



to create a convolution kernel K{r) that could be used to match the 
PSF of data with a width cti to a PSF with a width of a^. However, 
the MIPS 70 ^m data and the PACS 70 and 160 ^Im data have 
strong Airy rings, and applying this procedure would leave rings 
around bright sources. To correct for thi s, we created convo lution 
kernels following the proced ure given bv lGordon et"ai] ( l2008i) (and 
also seeh 



K 



following the procedure given by 
s iBendo et al] l l2010al) ). We used 

'W{lu)F[PSF2] 



F[PSFi] 



(3) 



to create convolution kernels K to match PSFi to PSF2. F is a 
Fourier transform and W{ij) is a radial Hanning truncation func- 
tion given by 



W{uj) 



uj > ujo 



(4) 



that is used to suppress high-frequency spatial noise in the result- 
ing kernels by adjusting uo. For the MIPS 70 ^m PSF, we used 
radially- smoothed versions of the empirical PSFs originally pre- 
sented in I Young etallJioO^ . For the PACS PSFs, we used radially- 
smoothed versions of the sum of three PSFs made from Vesta data 
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(see lLudl2O10L for more information]^- For the otlier PSFs, we 
used Gaussian functions. 

For quantitative analyses, we rebinned tiie data into 36 arc- 
sec pixels and then measured surface brightesses or intensities in 
all pixels where we detected emission at the 5a level or higher in 
the 70-500 fim bands. We used this bin size because it is equiva- 
lent to the FWHM of the PSF in the convolved data, and so we will 
avoid having multiple pixels sampling the PSF of a single point-like 
source. The coordinates of the bins were set up so that the centre 
of each galaxy would fall within the centre of a pixel. This binning 
procedure has an advantage over either measuring surface bright- 
nesses in regions identified by eye or measuring surface bright- 
nesses in regions identified in a specific wave band in that it will 
be unbiased towards any specific type of region and therefore will 
include both emission from point-like sources and the diffuse ISM. 
An example of the resulting convolved and binned data is shown 
in Figure |4] In our analysis, we use measurements in 322 bins for 
M81, 158 bins for M83, and 179 bins for NGC 2403. Note that we 
did not use the central 36 arcsec square bin in M81 for analysis be- 
cause it includes emission from the AGN, although we do display 
this data point in our figures. 



4 ANALYSIS 

4.1 Colour temperature maps 

Figures I5I7I show the colour temperatures based on the surface 
brightness ratios produced by this analysis. The colour tempera- 
tures are calculated using 

( f,hc/X,kT „ I 



+«9:15 



(5) 



in which dust emission is treated as originating from a single com- 
ponent that is optically-thin in the far-infrared. The parameter j3 
is the dust emissivity coefficient that indicates how the emissiv- 
ity scales as a fun ction of wavelength. We use j3 = 2 (given by 
iLi&DraindboOlh for this calculation, as it has previously been 
shown to accurately describe the observed emissivity of dust. The 
colour temperature maps for each galaxy look notably different in 
comparison to the stellar and dust structures in Figures 11131 The 
emission in some bands, particularly the shorter wavelength ones, 
may originate from dust with a range of temperatures. Moreover, 
temperatures determined from the Rayleigh- Jeans side of the SED 
may be improperly constrained without using data from the Wien 
side of the peak. Hence, these colour temperatures should be used 
to aid in relating the colour variations in the maps to the SED shape 
and should not necessarily be interpreted as the physical dust tem- 
peratures, even though this could be the case for the colour temper- 
atures based on the longer wavelength bands. Also, the 350/500 ^m 
colour temperatures tend to appear noisier. Both wave bands sam- 
ple locations on the Rayleigh-Jeans tail of the dust SED and are 
therefore relatively less sensitive to temperature variations, which 
causes the images to appear noisy. 

The M81 colour temperature map s are sim ply modified ver- 
sions of the ratio maps published bv ISendo ^et al. (2010b). The 
70/160 colour temperatures and, to a limited degree, the 



These data are available at ftp://ftp.sciops.esa.int/pub/hsc-calibration/ 
PACS/PSF/PACSPSF.PICC-ME-TN-033.vl.().tar.gz . The documentation 
is available from http://herschel.esac. esa.int/twiki/pub/Public/j 
PacsCalibrationWeb/bolopsfv 1.01 .pdf . 




+«8:50 



09:56 09:54 
Right Ascension (J2000) 



Figure 4. The original 250 /xm image of M81 (top), the 250 /xm image after 
it has been convolved with a kernel to match the PSF to that of the 500 /im 
data (middle), and the convolved 250 fjm image after rebinning (bottom). 
The green circles show the FWHM of the PSF of the data. In the convolved 
and rebinned maps, this PSF is 36 arcsec, which is also the size of the bins. 
Each map is 30 X 20 arcmin with north up and east to the left. 



160/250 /im colour temperatures are enhanced in the spiral arms. 
However, the 160/250 /im, 250/350 /im, and 350/500 /im colour 
temperature maps appear to be dominated by radial variations in 
the colour temperatures. This implies that the colour variations be- 
tween 160 and 500 /xm are dependent either on galactocentric ra- 
dius or on stellar surface brightness. 

The M83 colour temperature maps show enhancements near 
the nucleus and along spiral arms near star forming regions. The 
structures traced in the 70/160 /^m and 160/250 /im colour temper- 
ature maps clearly trace Ha regions along the edges of the spiral 
arms instead of the stellar structures seen in the 1.6 /im image, 
suggesting that these colour temperatures are linked to star forma- 
tion activity. We observe these structures on the leading edges of 
both spiral arms, which demonstrates that the structures are not 
simply the result of a misalignment in astrometry between two 
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Figure 5. The colour temperature maps of M81 and Ha and 1.6 fim images at the same resolution. All of these maps were created from data where the PSF 
was matched to the PSF of the 500 fim data, which has a FWHM of 36 arcsec (shown by the black circles in the lower left comer of each panel). This 36 arcsec 
circle is also the same width as the bins used in the analysis in Sections [4.2l and l43] Pixels where the signal was not detected at the 5o" level in either the single 
band shown (in the case of the Ha and 1.6 fim images) or in both bands (in the case of the colour temperature maps) were left blank (white). Each map is 
30 X 20 arcmin, and north is up and east is left in each panel. 



pairs of images. The relatively smooth gradient in the 250/350 nm 
colour temperature is consistent with the gradient in the 1.6 nm 
maps, and the central region with 250/350 fim colour temperatures 
above 30K, while very noisy, appears to correspond roughly to the 
bar seen in the 1.6 fim image. This suggests that the 250/350 fim 
colour temperatures could be affected by heating by the total stel- 
lar population. Foyle et al. (in preparation) are performing an anal- 
ysis that more carefully examines the colour gradients across the 
arms and the bar. While the nucleus itself appears to have enhanced 
70/160 fim colour temperatures, the regions with the most strongly 
enhanced colour temperatures in the other bands are the regions 
just to the east and west of the nucleus. The reason for this is un- 
clear. One possibility is that the PSF matching is imperfect and 
residual side lobes have appeared around the very bright nucleus 
in some of the data. However, such artefacts are not seen for other 



sources processed in the same way; see, for example, VS 44 in 
NGC 2403. Another possibility is that some of the dust in the nu- 
cleus itself is relatively cold because the dust shields itself from 
the nuclear starburst, but the regions immediately outside the nu- 
cleus are not as well shielded and thus become warmer in all wave 
bands. The nuclear colour temperature structure is poorly resolved 
in the resolutions that we are working at and the structures are rel- 
atively sensitive to 1-2 arcsec adjustments in astrometry (although 
the binned data are unaffected by such adjustments), so it is dif- 
ficult to accurately assess this with the data in hand. Future work 
with higher-resolution data is warranted. 

The NGC 2403 colour temperature maps show significant 
differences in the dust heating between 70 and 500 fim. The 
70/160 fim colour temperature map shows temperature enhance- 
ments mainly at locations with bright star forming regions. How- 
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Figure 6. The colour temperature maps of M83 and Ha and 1.6 ^m images at the same resolution. Each map is 20 X 20 arcmin. See the caption of Figure|5] 
for other information about the layout. 



ever, these regions become much less prominent in the colour tem- 
peratures based on longer wavelengths. In the 160/250 /im and 
250/350 /im colour temperature maps, the centre of the galaxy is 
approximately as warm as VS 44 (the bright star forming region to 
the northeast of the centre), and smooth radial variations generally 
become more prominent than localised dust temperature enhance- 
ments near star forming regions. The 350/500 /im map appears to 
be similar, although the map is so noisy that it is difficult to dis- 
cern any structures. Overall, these maps suggest that heating by star 
forming regions is important between 70 and 350 ^m but that vari- 
ations dependent on either total stellar surface brightness or radius 
becomes more important when proceeding to wavelengths longer 
than 160 fim, even though the dust structures themselves appear 
associated with the structures seen in Ha. 

Overall, the three galaxies show similar trends. In the colour 
temperatures based on the shortest wavelength data, the temper- 
atures appear strongly enhanced in regions with star forming re- 
gions, but in colour temperatures based solely on the >250 fim 
data, the variations in the colour temperatures appear more strongly 
related to either total stellar surface brightness or radius. This sug- 
gests a scenario in which star forming regions are located within 
dust clouds that are optically thick to the blue and ultraviolet light 
from star forming regions. In this scenario, the star forming regions 
would only heat a small fraction of the total dust mass in these 
clouds, and this warmer dust would be the predominant emission 
source at < 160 fim. However, most of the dust on the outside 



of these clumps as well as the dust in the diffuse ISM would be 
shielded from the ultraviolet and blue light from star forming re- 
gions and would instead be heated by the diffuse interstellar ra- 
diation field. While this dust would be cooler, the greater mass 
fraction of this cooler component would make it the more domi- 
nant emission source on the Rayleigh- Jeans side of the dust SED at 
> 250 fim. 

These imaging data have shown significant qualitatitve differ- 
ences between the dust heating as a function of wavelength. In the 
following sections, we will examine the different trends further us- 
ing quantitative analyses based on the binned data. 

4.2 Comparisons of surface brightness ratios to 1.6 /im 
emission 

Figure [8] shows how the surface brightness ratios for 36 arcsec 
square subregions vary versus 1.6 fim surface brightness, which ef- 
fectively traces the starlight from the total stellar populations within 
these galaxies. Statistical information on the relations is given in 
Table m 

Interestingly, strong correlations can be seen between the 
1.6 ^m surface brightness and the 160/250 fim, 250/350 ^m, and 
350/500 ^m ratios for all three galaxies. In almost all cases, the cor- 
relation coefficients are above 0.70. The squares of the correlation 
coefficients indicate the variance in one quantity that can be re- 
lated to the best fitting linear function to the other quantity. In this 
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Figure 7. The colour temperature maps of NGC 2403 and Hq and 1.6 fim images at the same resolution. Each map is 21 X 18 arcmin. See the caption of 
Figure|5]for other information about the layout. 



Table 4. Results from fits between surface brightness ratios and 1.6 /tm 
surface brightness. 



Galaxy 


Surface 


Slope" 


Correlation 




Brightness Ratio 




Coefficient" 


M81 


70/160 fim 


0.241 ± 0.002 


0.65 




160/250 tim 


0.142 ±0.004 


0.87 




250/350 Aim 


0.088 ± 0.005 


0.92 




350/500 fim 


0.067 ±0.007 


0.81 


M83 


70/160 


0.3456 ±0.0011 


0.82 




160/250 /xm 


0.1137 ±0.0005 


0.79 




250/350 fim 


0.0362 ± 0.0006 


0.83 




350/500 fim 


0.0474 ±0.0011 


0.64 


NGC 2403 


70/160 fim 


0.023 ± 0.002 


0.32 




160/250 /xm 


0.120 ±0.002 


0.85 




250/350 Aim 


0.074 ± 0.002 


0.92 




350/500 Aim 


0.047 ±0.003 


0.72 



" These quantities are for the relations describing the logarithm of the sur- 
face brightness ratios as a function of the logarithm of the 1.6 Atm surface 
brightness in Jy arcsec"^ . The fits are weighted by the uncertainties in the 
X- and y- values. 

context, at least 50% and often above 80% of the variance in the 
surface brightness ratios can be related to the stellar surface bright- 
ness. The 350/500 A^ni ratios for the three galaxies show more scat- 
ter (i.e. have lower correlation coefficients) than the 250/350 A^m 
ratios in the relations with 1.6 fim surface brightness. This could 



be because of the relatively high uncertainties in the ratios, which 
are in part a consequence of sampling data on the Rayleigh- Jeans 
tail of the SED where temperature variations have a relatively small 
effect on the colours. Additionally, the nucleus of M81 (which cor- 
responds to the data point with the highest 1.6 Aim surface bright- 
ness) deviates significantly from the relations between the 1.6 
surface brightnes s and either the 250 /350 A^rn or 350/500 A^rn ra- 
tios. As stated bv lBendoetai] ( l2010bl) . this is probably because the 
AGN produces non-t hermal emission at submillimetre wavelengths 
jMarkoff et al.ll20o3) . and it was not included in the derivations of 
the statistical results shown in Table [8] None the less, the over- 
all results indicate that the total stellar populations in each galaxy 
play a significant role in locally heating the dust that is emitting 
at >160 Atm. Such a result may be expected in M81, where the 
bulge has a significantly higher surface brightness in visible light 
compared to the disc and where dust heating could be expected to 
be dominated by evolved stars. However, it is surprising that this 
strong correlation is seen in NGC 2403, where the bulge is either 
very faint or non-existent. 

The results for the correlations between the 70/160 A^m ra- 
tios and the 1.6 fim surface brightness, however, show that dust 
heating by the total stellar population may not be as important, al- 
though this varies among the three galaxies examined here. Weak 
correlations are seen in M81; less than 50% of the variance in the 
ratio is related to the stellar surface brightness. A relation may be 
present where the 1.6 A^m surface brightness is higher than 10~* 
Jy arcsec^ in M81, but the relation flattens out for fainter regions. 
This could indicate that the bulge plays a role in heating the dust 
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Figure 8. Comparisons of the surface brightness ratios to the 1.6 fim surface brightnesses for the 36 arcsec square subregions in the sample galaxies. The solid 
lines show the best fit lines. The thick red line shows the weighted mean of the ratios averaged over intervals with a width of 0.2 in log(7i/(1.6 fim)), and the 
thinner red lines show the weighted standard deviation of the data (except when one data point falls within an interval, in which case the uncertainty in the 
data point is used). Parameters related to the fits are shown in Table|4] 



that produces the 70 emission, but only in the centre where 
the bulge is very bright. The correlation coefficient for the relations 
between the 70/160 fim ratio and Ha (discussed in Section l43t 
is higher for M83, and the 70/160 /im colour temperature map 
for M83 shows that the regions with enhanced temperatures cor- 
respond to Ha sources, which suggests that the relation seen be- 
tween the 70/160 fj.m ratio and 1.6 /xm surface brightness is partly 
a consequence of the correlation between the stellar structures and 
the structures tracing star formation. As for NGC 2403, the corre- 
lation between the 70/160 ^m ratio and 1.6 /im surface brightness 
is negligible, which reflects the weak influence of the total stellar 
population on the colour temperatures, as we also inferred from 
Figure [T] 



iBendo et al. I( l2010bh noted that the far-infrared surface bright- 
ness ratios for M81 were correlated with radius. This result was 
used to infer that the ratios depended on the total stellar surface 
brightness, which varies with radius within nearby galaxies. For 
completeness, we re-examined the relations for all three galaxies in 
this study. Figure |9] shows how the surface brightness ratios for the 
36 arcsec square subregions vary versus deprojected galactocentric 
radius, and Table[5]gives statistical information on these relations. 
Visually, the relations between radius and either the 160/250 /im, 
250/350 fim, or 350/500 ^m ratios appear almost as strong as the 
relations between 1.6 ^im surface brightness and these ratios. This 
is mainly because the stellar surface brightness itself is correlated 
with radius. However, the correlation coefficients of the trends for 
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Table 5. Results from fits between surface brightness ratios and deprojected 
galactocentric radius. 



Galaxy 


Surface 


Slope" 


Correlation 




Brightness Ratio 




Coefficient" 


M81 


70/160 ^im 


-0.0480 ± 0.0005 


-0.57 




160/250 Mm 


-0.0268 ± 0.0007 


-0.82 




250/350 iim 


—0 01 fi4 + oooq 


-0.92 




350/500 Aim 


-0.0121 ±0.0013 


-0.80 


M83 


70/160 fim 


-0.1122 ±0.0003 


-0.77 




160/250 Atm 


-0.0331 ±0.0001 


-0.76 




250/350 Aim 


-0.0112 ±0.0002 


-0.75 




350/500 Aim 


-0.0142 ± 0.0003 


-0.66 


NGC 2403 


70/160 Aim 


-0.0084 ± 0.0007 


-0.26 




160/250 Atm 


-0.0350 ±0.0006 


-0.89 




250/350 Atm 


-0.0221 ±0.0006 


-0.89 




350/500 Aim 


-0.0157 ±0.0010 


-0.80 



" These quantities are for the best fit line describing the logarithm of the 
surface brightness ratios as a function of deprojected galactocentric distance 
in kpc. The slope has units of dex kpc~^. 



the relation with 1.6 Aim surface brightness have slightly higher ab- 
solute values for all of the M81 relations, three of the relations for 
M83, and two of the NGC 2403 relations. Moreover, the radial pro- 
files for M81 depart from the best fit lines in the central 2 kpc in 
a way that would be consistent with heating by bulge stars in this 
region. Furthermore, the colour temperature maps in Figures ISTTI 
show the presence of some non-axisymmetric structures that cor- 
respond to spiral arms or filaments in each of the galaxies. Some 
of these features correspond to locations with star forming regions, 
which can be thought of as also enhancing the total stellar surface 
brightness. Other features correspond to structures in the 1.6 Atm 
image that are more clearly cases where the surface brightness of 
the total stellar population increases. Hence, it would make more 
sense to infer that the surface brightness ratios are more strongly 
dependent on stellar surface brightness than radius. 

Physically, a relation between stellar surface brightness and 
dust colours is more straightforward than a relation between galac- 
tocentric radius and dust colours. The stars can directly heat the 
dust, so the stellar surface brightness would be directly linked to 
dust temperature (and hence the surface brightness ratios) regard- 
less of the relation between stellar surface brightness and radius. 
However, to link the dust temperature with radius without connect- 
ing dust temperature to stellar surface brightness is very difficult. 
We therefore conclude that the correlations between galactocentric 
radius and either the 160/250 Aim, 250/350 Aim, or 350/500 Aim ra- 
tio actually reflect radial variations in the starlight heating the dust. 

Figure[TO]shows how the 70-500 Aim surface brightnesses for 
the 36 arcsec square subregions are relate to the 1.6 A*m surface 
brightnesses. This comparison shows how the dust emission may be 
related to the distribution of the total stellar population. Good cor- 
relations, with coefficients above 0.89, are found between the dust 
and stellar emission for M83, which we expected, as the infrared 
and stellar emission appear qualitatively similar. Surprisingly good 
correlations, with correlation coefficients above 0.85, are found be- 
tween dust and stellar emission in NGC 2403 despite the signficant 
diffrences in the appearance of the stellar and dust emission. How- 
ever, both the stellar and dust emission generally decreases with 
galactocentric radius, which could cause the stellar emission to ap- 
pear correlated overall with the dust emission. In contrast to the 
other two galaxies, the relations between stellar and dust emission 



Table 6. Results from fits between far-infrared and 1.6 A*m surface bright- 
nesses. 



Galaxy 


Infrared Surface 
Brightness 


Slope" 


Correlation 
Coefficient" 


M81'' 


70 A^m 


0.4164 ± 0.0015 


0.78 




160 Aim 


0.1946 ± 0.0019 


0.69 




250 Aim 


0.115 ± 0.003 


0.52 




350 Aim 


0.077 ± 0.004 


0.34 




500 Aim 


0.048 ± 0.006 


0.15 


M83 


70 Aim 


1.19 ± 0.05 


0.90 




160 Aim 


0.84 ± 0.03 


0.91 




250 Aim 


0.9284 ± 0.0015 


0.90 




350 Aim 


0.8745 ± 0.0014 


0.90 




500 Aim 


0.7925 ± 0.0018 


0.89 


NGC 2403 


70 Aim 


0.8038 ± 0.0016 


0.88 




160 Aim 


0.649 ± 0.003 


0.90 




250 Aim 


0.5476 ± 0.0018 


0.89 




350 Aim 


0.462 ± 0.002 


0.87 




500 Aim 


0.389 ± 0.003 


0.87 



" These quantities are for the relations for the logarithms of the surface 
brightnesses in Jy arcsec^^. 



in M81 is fairly good at 70 and 160 Aim but worsens towards longer 
wavelengths. Emission in the 70 and 160 Aim bands is more sensi- 
tive to variations in dust temperature, and since these bands partly 
trace dust within M81 that appears to be heated by the stars seen at 
1.6 Aim, the 70 and 160 Aim emission is correlated with the 1.6 Aim 
emission. However, variations in the 350 or 500 Aim emission may 
be more dependent on dust surface density than dust temperature, 
and so these wave bands appear more poorly correlated with the 
stars in M81 because the stellar mass is distributed differently from 
the dust mass. 

Although we have found that the distribution of dust and stel- 
lar emission may match in many circumstances, this alone does 
not necessarily indicate that the dust traced in these wave bands 
is heated by the total stellar population. Dust emission in a single 
wave band is a function of both dust temperature and dust surface 
density. For example, it is possible for dust emission at 70 Aim in 
NGC 2403 to appear correlated with 1.6 Aim emission even though 
the colour temperature map and the analysis of the surface bright- 
ness ratios suggest that the dust is not heated by the total stellar 
population. Instead, the dust and stellar surface densities are both 
correlated with radius in this galaxy, causing the 70 Aim emission 
to appear related to 1.6 Aim emission. Conversely, it is possible for 
dust emission in a single wave band, such as the 500 Aim data for 
M8 1 , to appear uncorrelated with stellar surface brightness because 
the wave band is more strongly affected by dust surface density 
variations than dust temperature variations and because the dust 
and stellar mass is distributed differently within the galaxy. None 
the less, this does not preclude the possibility that the dust is heated 
primarily by the total stellar population, thus causing the tempera- 
ture of the dust to be dependent upon the stellar surface brightness. 



4.3 Comparisons of colour temperatures to star formation 

Figure [TT] shows how the colour temperatures for 36 arcsec square 
subregions vary versus Ha intensity, which traces the star forma- 
tion within the galaxies. The statistics from these relations are given 
in Table|7] Each of the galaxies in this analysis give intriguingly dif- 
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Figure 9. Compaiisons of the surface brightness ratios to deprojected galactocentric radius for the 36 arcsec square subregions in the sample galaxies. The 
solid lines show the best fit lines. The thick red line shows the weighted mean of the ratios averaged over intervals with a width of 1 kpc, and the thinner red 
lines show the weighted standard deviation of the data (except when one data point falls within an interval, in which case the uncertainty in the data point is 
used). Parameters related to the fits are shown in Table [5] 



ferent results, especially when the results are compared with those 
data shown in Figure[8]and Table|4] 

The data for M81 show that the surface brightness ratios are 
generally poorly correlated with Ho intensity. This is also borne out 
by the general absence of spiral structure in the colour temperature 
maps in Figure [5j the colour temperatures primarily vary with ra- 
dius, and, with the exception of Tc (70/160 ^m), the spiral structure 
is very weakly enhanced in the spiral arms if it is enhanced at all. 
It is possible, however, that the 70/160 ^m surface brightness ratio 
is enhanced in star forming regions in the spiral arms, as the re- 
gions appear relatively bright in the colour temperature map and as 
some of the data with the highest Ha values in Figure[TT]do appear 
to show a relation with the 70/160 /im ratio. However, the central 



~ 3 kpc of M81, which is devoid of star formation and which has 
low Ha intensities (although incompletely-subtracted continuum 
emission from the bulge is seen in the Ha image) exhibits higher 
70/160 ^m ratios, demonstrating that the evolved stellar population 
must be the heating source for the dust emitting at 70 /im in this 
region. As for the other surface brightness ratios, the stronger cor- 
relations between the ratios and 1.6 ^m surface brightness clearly 
show that the total stellar population plays a much larger role than 
the star forming regions in heating the dust in M8 1 . 

For M83, the correlation coefficients for the relations between 
Ha intensity and either the 70/160 ^im or 160/250 /im ratios is 
slightly higher than for the corresponding relations with 1.6 ^m 
surface brightness, but the correlation coefficient for the Ha in- 
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Figure 10. Comparisons of infrared surface brightnesses in individual wave bands to the 1 .6 fim surface brightnesses for the 36 arcsec square subregions in 
the sample galaxies. The solid lines show the best fit lines. The thick red line shows the weighted mean of the ratios averaged over intervals with a width of 0.2 
in log(7y(1.6 /im)), and the thinner red lines show the weighted standard deviation of the data (except when one data point falls within an interval, in which 
case the uncertainty in the data point is used). Parameters related to the fits are shown in Table|6] 



tensity and either the 250/350 or 350/500 fim ratios is lower 
than the corresponding relations with 1.6 /im surface brightness. 
For the 70/160 fim, 160/250 ^m, and 250/350 fj,m ratios, the co- 
efficients for the relations with Ha intensity are within 0.10 of 
the coefficients for the relations with 1.6 fim surface brightness. 
This implies that both 1.6 ^tm and Ha emission correlate with dust 
temperatures almost equally well. This could be true for a couple 
of reasons. First of all, the Ha and 1.6 /im emission both trace 
very similar spiral structures within M83 and are correlated with 
each other to some degree, especially when the data are binned 
up as we have done here. Hence, even if Ha or 1.6 fim emission 
trace completely different stellar populations, a relation found be- 
tween either band and far-infrared colours may naturally apply to 



the other band as well. It is also possible that both star forming re- 
gions and evolved stars heat the dust observed in the far-infrared, 
with evolved stars responsible for a slightly greater fraction of the 
heating for dust emitting at longer wavelengths and star forming 
regions heating more dust at shorter wavelengths. A third possi- 
bility is that the 1.6 fim band includes a significant contribution 
from young star-forming regions. However, the spiral structure in 
the 1.6 fim image is not as strong as the structure in the Ha im- 
age, and both the 1.6 fim and Tc (250/350 fim) images exhibit bar 
structures that do not have such apparent counterparts in either the 
Ha, Tc (70/160 fim), or rc( 160/250 fim) images, implying that 
"contamination" of the 1.6 fim band by young stars cannot be the 
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Figure 11. Comparisons of the surface brightness ratios to the Ha intensities for the 36 arcsec square subregions in the sample galaxies. The solid lines show 
the best fit lines. The thick red line shows the weighted mean of the ratios averaged over intervals with a width of 0.2 in log(7(Ho)), and the thinner red 
lines show the weighted standard deviation of the data (except when one data point falls within an interval, in which case the uncertainty in the data point 
is used). Data from the inner central 5.7 X 3.0 arcmin (approximately the central 3 kpc) of M81 are not included in the analysis, as it primarily includes 
incompletely-subtracted continuum emission from the bulge. Parameters related to the fits are shown in Table|7] 



only explanation for the similarities in the correlations of surface 
brightness ratios with either 1.6 or Ha emission. 

The results for NGC 2403 are similar to M83 in that the 
70/160 /im are clearly more strongly correlated with Hq inten- 
sity but the 250/350 fim surface brightness ratio is more strongly 
correlated with 1.6 fim surface brightness, with the correlation co- 
efficients for the relations between the 160/250 /im ratio and ei- 
ther Ha intensity or 1.6 /im surface brightness appearing almost 
equal. (The correlations between the 350/500 /im ratio and either 
Ha intensity or 1.6 fim surface brightness appear almost equally 
correlated probably because of the noisiness of the ratio.) The dif- 
ference in the correlation coefficients for the relations between the 
250/350 /im ratio and either 1.6 fim or Ha emission would not by 



itself necessarily indicate that the total stellar population is a sig- 
nificantly stronger source of heating the dust at these wavelengths. 
However, the similarities in these correlation coefficients may be 
the consequence of the high scatter at low surface brightness levels 
in the data. An examination of the regions with the highest surface 
brightnesses show that the 250/350 fim ratios are more strongly re- 
lated to the 1.6 nm surface brightnesses. The 250/350 fim ratios 
for the binned data in Figure [8] fall within 2a of the best fitting 
relation with the 1.6 /im surface brightness for the regions where 
1^(1.6 fj.m> 10^'' Jy arcsec^^. In contrast, the 250/350 /im ratios 
for data with the highest Ha intensities (where/ (Ha > 10^^^ erg 
cm~^ s~^ arcsec"^) in Fi gure [TTI frequently fall > 3a away from 
the best fit line. This variation in the dispersion can also be seen 
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Table 7. Results from fits between surface brightness ratios and Ha inten- 
sities. 



Table 8. Results from fits between infrared surface brightnesses and Ha 
intensities. 



Galaxy 


Surface 
Brightness Ratio 


Slope" 


CoiTelation 
Coefficient" 


Galaxy 


Infrared Surface 
Brightness 


Slope" 


Correlation 
Coefficient" 


M81'' 


70/160 


0.178 ± 0.005 


0.52 


M81'' 


70 ^m 


0.834 ± 0.005 


0.80 




160/250 /jm 


0.088 ± 0.006 


0.45 




160 fim 


0.660 ± 0.005 


0.74 




250/350 /im 


0.037 ± 0.007 


0.33 




250 fivn 


0.523 ± 0.006 


0.78 




350/500 /im 


0.039 ± 0.010 


0.33 




350 fj.m 


0.485 ± 0.006 


0.79 


M83 


70/160 fim 


0.3219 ± 0.0006 


0.89 




500 /im 


0.433 ± 0.009 


0.78 




160/250 fj.m 


0.1174 ±0.0005 


0.87 


M83 


70 fim 


0.95 ± 0.03 


0.93 




250/350 fim 


0.0352 ± 0.0006 


0.77 




160 fim 


0.66 ±0.02 


0.92 




350/500 ^m 


0.0425 ± 0.0011 


0.50 




250 /im 


0.8018 ± 0.0006 


0.90 


NGC 2403 


70/160 ^im 


0.144 ±0.002 


0.51 




350 /Jim 


0.7608 ± 0.0017 


0.89 




160/250 ^lm 


0.147 ±0.003 


0.86 




500 /im 


0.7153 ±0.0012 


0.90 




250/350 fim 


0.064 ± 0.002 


0.78 


NGC 2403 


70 fim 


0.7474 ± 0.0012 


0.91 




350/500 ^tm 


0.051 ± 0.003 


0.74 




160 fim 


0.600 ± 0.003 


0.88 












250 fim 


0.511 ± 0.002 


0.86 


" These quantities are for the relations describing the lof 


>arithm of the sur- 




350 fim 


0.447 ± 0.002 


0.86 


face brightness ratios as a function of the logarithm of the Ho intensity in 

_9 _1 _9 




500 fim 


0.394 ± 0.003 


0.86 



erg cm s arcsec . 

^ Data from the inner central 5.7 X 3.0 arcmin (approximately the cen- 
tral 3 kpc) of M81 are not included in the analysis, as it primarily includes 
incompletely-subtracted continuum emission from the bulge. 




0.010 0.100 1.000 
l(normalized) 



Figure 12. The weighted standard deviation of the 250/350 fim surface 
brightness ratio for the 36 arcsec square subregions in NGC 2403 plotted as 
a function of the Ha intensity (blue) and 1.6 fim surface brightness (red). 
These curves are the same values used to plot the weighted standard devia- 
tions in Figures [8] and [TT] The Ha intensity and 1.6 fim surface brightness 
have been normalised by the value of the highest surface brightness interval. 



" These quantities are for the relations describing the logarithm of the sur- 
face brightnesses in Jy arcsec as a function of the logarithm of the Ha 
intensity in erg cm^^ s~^ ai'csec"^. 

Data from the inner central 5.7 X 3.0 arcmin (approximately the cen- 
tral 3 kpc) of M8 1 ai'e not included in the analysis, as it primarily includes 
incompletely-subtracted continuum emission from the bulge. 



We also used our binned data to compare surface brightnesses 
from individual far-infrared bands to Ha emission so as to attempt 
to replicate the relations between star form ation tracers and 100 , 
160 , or 250 fim emission in M33 reported bv lBoguien et alj ( l201(]|) 
and lVerlet et"aD (2010). Fi gure[T3] shows the result of this compar- 
ison, and Table [8] shows the slopes and correlation coefficients for 
these relations. The correlations are notably strong. In all cases, 
the correlation coefficients are above 0.7, which indicates that at 
least > 50% of the variance in the emission in any of these infrared 
bands can be accounted for by the relation with Ha emission. These 
strong correlations are even seen in the relations between Ha inten- 
sity and 500 fim surface brightness. However, paired with the re- 
sults on how the surface brightness ratios vary with 1.6 fim surface 
brightness and Ha intensity, it is clear that the relations seen in Fig- 
ure[T3]are not necessarily a consequence of the dust being directly 
heated by the star forming regions. Instead, the relations shown in 
Figure[T3]may arise as a result of the dust being located in the same 
regions as the star forming regions. We discuss the further implica- 
tions of the relations between infrared surface brightness and Ha 
intensities in the next section. 



in Figure [12] which shows the dispersion in the 250/350 fim ratio 
as a function of either the Ha or 1.6 fim emission (after normal- 
ising both to the peak surface brightness). Additionally, the quali- 
tative similarities between the Tc (250/350 fim) and 1.6 fim maps 
and the notable differences between the Tc (250/350 fim) and Ha 
maps in Figure [7] indicate that the total stellar population is more 
important than star forming regions for heating the dust observed 
at > 250 fim. While such a result may not be surprising for M81, 
where the bulge is very large compared to the disc, this is particu- 
larly surprising for NGC 2403, which is a relatively bulgeless late- 
type spiral galaxy with widespread star formation in its disc, and 
we anticipate that these results may apply to other late-type spiral 
galaxies as well. 



4.4 Relative contributions of total stellar populations and 
star formation to dust heating 

In prior sections, we performed comparisons between surface 
brightness ratios and the emission from either the total stellar pop- 
ulations or star forming regions. When the surface brightness ra- 
tio variations are much more strongly correlated with one source 
of dust heating than with the other, we can draw very sharp con- 
clusions about the dust heating. For example, we can confidently 
conclude that the 70/160 fim surface brightness ratio variations in 
NGC 2403 are primarily dependent on heating by star forming re- 
gions, while we can also say that the 250/350 fim variations in 
M81 are dependent on heating by the total stellar population. How- 
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Figure 13. Compaiisons of infrared surface brightnesses in individual wave bands to tlie Ha intensities for the 36 arcsec squai'e subregions in the sample 
galaxies. The solid lines show the best fit lines. The thick red line shows the weighted mean of the ratios averaged over intervals with a width of 0.2 in 
log(/(Ha)), and the thinner red lines show the weighted standard deviation of the data (except when one data point falls within an interval, in which case the 
uncertainty in the data point is used). Data from the inner central 5.7 X 3.0 arcmin (approximately the central 3 kpc) of M81 are not included in the analysis, 
as it primarily includes incompletely-subtracted continuum emission from the bulge. Parameters related to the fits are shown in Table[8] 



ever, the above comparisons have not been as effective in revealing 
which dust heating source is dominant in cases where the corre- 
lations produce similar correlation coefficients, as, for example, in 
the case of the 160/250 ^m data for NGC 2403. 

To address this issue, we attempt to fit the surface brightness 
ratios for the 36 arcsec square subregions as a function of both 
1.6 fim and Ha emission. The equation we fit to the data is based 
on 

CT^+^^Eir,, (6) 

a generalised version of the Stefan-Boltzmann law for the case of 
blackbody emission modified by an emissivity function described 
by a power law with a coefficient /3. Ejn will be equal to the total 



energy from all dust heating sources and can be written as the sum 
of scaled versions of the 1.6 fj,m and Ha emission. For a small 
range in temperature variation (as is the case for our data), T can 
be related to a surface brightness ratio through a power law. The 
difference between the actual colour temperatures and the colour 
temperatures from this power law differ by only ~ 2% over the 
range of surface brightness ratios that we are using; this is ususally 
lower than the uncertainties in the ratios themselves. Given this, we 
can represent Equation|6]as 

Ci 1 t4t4 ) ° = C2/(Ha) + CsLilMfim), (7) 
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or we can rewrite this in logaritlimic form as 

In = «ln(C2/(Ha) + C3h{1.6^im)) - In(Ci). (8) 

Equation[8]can be simplified as 

ln( :^4t4 ) =«ln(/(HQ)+Ai/,(1.6/im)) + A2. (9) 

Equation |9] can be thought of as similar to the relations typically 
shown between colour and flux or variants of these quantities (such 
as Figures Island I lit but with two flux quantities simultaneously 
fit to the colour instead of one. The relative contribution of star 
forming regions to the dust heating can be calculated using 



Erotai /(Ha) +Ai7,(1.6Aim)' 

This approach to determining the relative contributions of star 
forming regions and the total stellar population to dust heating is 
relatively simplistic. More complex models would be able to more 
accurately describe the propagation of photons through the ISM, 
the absorption of light by dust, the range of dust temperatures, 
and the dust emission, but this requires making many assumptions 
about the stellar populations and the physical properties of the dust 
particles. In contrast, the analysis we have outlined above relies 
primarily upon the assumptions that dust heating occurs locally and 
that the relation between surface brightness ratios and temperatures 
follow power laws. Hence, it can be used as a way to examine dust 
heating from two energy sources using measurements that are inde- 
pendent of dust surface density in a way that relies upon very few 
assumptions. 

We fit Equation [8] to the 36 arcsec square subregions in each 
galaxy, as it was easier to use this equation with a Levenberg- 
Marquardt algorithm to converge upon the best fit to the data, even 
though the solution is degenerate. We then used these results to 
calculate the parameters for Equation |9] Table |9] given the best fit- 
ting parameters along with the correlation coefficients for the val- 
ues given by the right and left sides of Equation |9] and the frac- 
tional contribution of star forming regions to dust heating as given 
by Equation (To) Uncertainties in the quantities in Table [9] are cal- 
culated using a monte carlo approach in which we perform fits in 
a series of iterations by adding Gaussian noise to the data and then 
measure the standard deviation in the results. Figure [14] shows the 
relations between the surface brightness ratios and the best fitting 
relations. Additionally, we used these results along with Equation|5] 
to create colour temperature maps constructed from the Ha and 
1.6 fj-m images. These reconstructed colour temperature maps as 
well as the colour temperature maps based on the observed surface 
brightness ratios and maps of the relative contribution of star form- 
ing regions to dust heating are all shown in Figures [l5ll7l 

We encountered two problems with fitting Equation [8] to the 
data. In fitting the 70/160 /im surface brightness data for M83, we 
found that the difference between the ratio and the resulting fit was 
significantly high around the nucleus, and the reconstructed colour 
temperature map in Figure [T6]did not look like the observed colour 
temperature map. Additionally, the fractional contribution of star 
formation to dust heating seemed abnormally low; when the nu- 
cleus is included, the fractional contribution of star formation to 
dust heating is ~ 0.3. This implies that the nucleus does not fol- 
low the same 70/160 fim dust heating relation as the rest of the 
disc, for reasons that are unclear but that could be related to hot 
(~ 100 K) dust heated by the nuclear starburst that is visible in 
the 70 /im band. Therefore, we excluded the central 3x3 regions 



(108 X 108 ai-csec) of M83 for the analysis on the 70/160 fim ratio. 
(Excluding the same central region when fitting Equation [8] to the 
160/250, 250/350, or 350/500 /im data did not affect the fits signfi- 
canfly.) In fitting Equation[8]to the 70/160 /im data for NGC 2403, 
we found that the 6*3/^(1. 6/im) term was negligible compared to 
the C2/(Ha) term, which caused convergence problems when try- 
ing to perform a nonlinear least squares fit. We therefore report 
results in which the 70/160 /im ratio is a function of only the Ha 
intensities (which is consistent with the interpretation that varia- 
tions in dust heating traced by this ratio depend primarily on star 
formation). 

After dealing with these two exceptions, the resulting fits are 
generally very good. The correlation coefficients in Table[9]for the 
relations between the logarithm of the surface brightness ratios and 
the expression on the right side of Equation |9] are generally higher 
than the corresponding correlation coefficients in either Tables|4]or 
[7] for fits to either the 1.6 fim or Ha emission, although in many 
cases, the improvement in the correlation coefficients is marginal. 
Moreover, most correlation coefficients in Table |9] are above 0.7, 
indicating that at least 50% of the variance in the data can be ac- 
counted for by the relations. The only exceptions are the 70/160 /im 
ratio for NGC 2403, which we identified as problematic up above, 
and the 350/500 /im ratio for M83, which is very noisy in general. 

The reconstructed colour temperature maps based on Equa- 
tion[9]in Figures [l5II7l look qualitatively similar to the colour tem- 
perature maps based on the observed surface brightness ratios. The 
reconstructed maps generally match the overall radial gradients in 
temperature seen in many of the observed maps, and the recon- 
structed maps also show that we could reproduce the colour tem- 
peratures of many of the individual star forming regions within the 
spiral structures of each galaxy. However, the temperatures of the 
reconstructed star forming regions do not always precisely match 
the colour temperatures in the observed map. This could be be- 
cause of variable extinction among individual star forming regions, 
which could cause the relation between star formation (as traced 
by Ha emission) and dust heating to vary. We also occasionally see 
some residual structures in temperature that show features that are 
not quite described by the reconstructed map. In M8 1 , the nucleus 
(excluded from the fits) appears notably cold in the residual maps 
because of the synchrotron emission from the AGN. The nucleus 
of M83, which caused problems in the fits and was excluded, ap- 
pears as a hot spot in the residual 70/160 fim map. Residual spiral 
structure is seen in the 160/250 /im and 250/350 /im residual maps, 
which could be indicative of dust heating gradients across the spi- 
ral arms that are inadequately described by Equation |9] and show 
the limitations of using this approach to reproducing dust colour 
temperature variations. Foyle et al. (in preparation) found a simi- 
lar offset between dust temperatures and star formation in the spi- 
ral arms. Finally, we see residual structure in the 70/160 /im map 
of NGC 2403 that is clearly a result of latent image effects in the 
70 /im data. 

The main purpose of these fits was to determine the relative 
contributions of star forming regions and the total stellar popula- 
tions to dust heating. The results in Table |9] generally show a de- 
crease in the relative contribution of dust heating from star forma- 
tion as wavelength increases. In M83 and NGC 2403, star forming 
regions contribute > 50% of the total heating of the dust traced 
by the < 250 /im bands, although in M81, the dust heating up to 
70 fim appears to still be dominated by the total stellar population. 
At longer wavelengths, the relative contribution of star forming re- 
gions decreases signficantly. Star forming regions appear to con- 
tribute only < 30% of the energy for heating the dust traced by 
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Figure 14. Comparisons of the surface brightness ratios to the relation described by Equation |2]for the 36 arcsec square subregions in the sample galaxies. 
The solid lines show the best fit lines, which have a slope of 1. The thick red line shows the weighted mean of the ratios averaged over different interval 
widths for each surface brightness ratio (0.2 in the units of the x-axes for the 70/160 ^m ratio, 0.1 for the 160 /im ratio, and 0.03 for the 250/350 fim and 
350/500 fim ratios), and the thinner red lines show the weighted standard deviation of the data (except when one data point falls within an interval, in which 
case the uncertainty in the data point is used). Data from the inner central 3x3 subregions (108 X 108 arcsec) of M83 are not included in the analysis, as it 
caused problems during the fit as described in the text. Parameters related to the fits are shown in Table [9] 
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Figure 15. Ttie left column shows the observed colour temperature maps of M8 1 . The next column shows the colour temperature maps reconstructed using 
Equation [s] and |9] the Ha and 1.6 fim images, and the parameters if Table|2] The third column shows the residual values when the reconstructed image is 
subtracted from the observed image. The fourth column shows the fraction of dust heating that is attributed to star formation as traced by the given surface 
brightness ratio; these images are created using Equation[To] the Ha and 1.6 fim images, and the Ai parameters in Table|2] All of these maps were created 
from data where the PSF was matched to the PSF of the 500 fim data, which has a FWHM of 36 arcsec (shown by the black circles in the lower left comer of 
each panel). This 36 arcsec circle is also the same width as the bins used in the analysis in Sections [4. 2l and l4!3] Pixels where the signal was not detected at the 
5<T level in both bands used to create the observed colour temperature maps were left blank (white). In the other three columns, these pixels are left blank and 
pixels where data are not detected at the 5a level in either the Ha or 1.6 fim images were left blank. Each map is 30 X 20 arcmin, and north is up and east is 
left in each panel. 
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Table 9. Results from fitting surface brightness ratios as a function of both Ha and 1 .6 fim emission. 





SurfflCG 
Brightness Rfltio 






^1 


^2 






C orre 1 sti on 


M81 


70/160 ^tm 


(2.25 - 


1-/ 


- 0.05) X 10-12 


12.01 ±0.12 


0.369 ± 0.004 


0.361 ±0.006 


0.76 




160/250 ^m 


(5.4- 


h/ 


- 0.6) X 10-12 


7.3 ±0.2 


0.185 ±0.005 


0.19 ±0.02 


0.89 




250/350 ^m 


(1.3- 


h/ 


- 0.4) X 10-" 


4.3 ±0.2 


0.100 ±0.006 


0.09 ± 0.03 


0.93 




350/500 ^m 


(1.2- 


VI 


- 1.1) X 10-" 


3.6 ±0.3 


0.078 ± 0.007 


0.10 ±0.06 


0.82 


M83 


70/160 ^im 


(5.7 4 


-1 


- 1.2) X 10-13'= 


8.91±0.12<= 


0.274 ±0.004= 


0.924 ± 0.014= 


0.89= 




160/250 Atm 


(1.06 - 


VI 


- 0.08) X 10-12 


5.21 ±0.02 


0.120 ±0.001 


0.854 ± 0.009 


0.88 




250/350 ^m 


(1.02 - 


VI 


- 0.17) X 10-11 


2.24 ±0.03 


0.038 ±0.001 


0.38 ±0.04 


0.85 




350/500 fim 


(1- 


VI 


-3) X 10-10 


2.63 ±0.06 


0.048 ± 0.001 


0.06 ±0.04 


0.66 


NGC 2403 


70/160 fj,m 






0/ 


4.23 ± 0.02-^ 


0.144 ±0.002-f 


1/ 


0.51-f 




160/250 ^lm 


(6.2 - 


VI 


- 0.6) X 10-12 


6.39 ±0.10 


0.158 ±0.003 


0.65 ±0.02 


0.89 




250/350 Atm 


(1.1- 


VI 


- 0.3) X 10-10 


3.32 ±0.09 


0.077 ± 0.002 


0.10 ±0.03 


0.93 




350/500 ^m 


(1.6- 


VI 


- 0.6) X 10-11 


2.83 ±0.13 


0.055 ±0.002 


0.42 ± 0.08 


0.76 



" These parameters correspond to the ones in Equation|9] 

The units of Ai are (erg cm-2 s— l arcsec-2) (MJy sr-i)— l. 
= This is relative fraction of dust heating from star forming regions as given by Equation llOl The value is calculated with the data for all 
36 arcsec square regions used when fitting Equation|9]to the data. 

These are correlation coefficients for the relations between the right and left sides of Equation|9] 
= These fits were performed excluding the central 3 X 3 36 arcsec square regions. See the text for details. 

f Instead of reporting the result from the fit to Equation[8]here, we instead report the result of the fit of a power law relation between the 
70/160 fim surface brightness ratio and the Ha intensity. See the text for details. 



the 250/350 fim and 350/500 fim bands in all three galaxies. The 
relation for the 350/500 fim ratio for NGC 2403 seems anomalous 
in that it implies that ~ 40% of the energy for the dust heating 
originates from star forming regions, which is higher than what is 
implied by the fit to the 250/350 ^m ratio for that galaxy. How- 
ever, we suspect that this is because the 350/500 ^m data are on the 
Rayleigh- Jeans side of the SED and are strongly affected by noise. 
Hence, colour variations related to dust heating are not detected at 
a high signal-to-noise level in NGC 2403, and so the methods we 
are using may not be effective in this one case. 

Overall, the results here are consistent with the interpreta- 
tions we have provided for the results in Sections 14.21 and 14.31 At 
wavelengths shorter than 160 /im, most of the dust appears to be 
heated by star forming regions in these galaxies, which is why the 
70/160 fj,m ratios appear more closely dependent on the Ha surface 
brightnesses. At wavelengths longer than 250 ^m, the dust appears 
to be primarily heated by the total stellar populations, and so the 
250/350 fim and 350/500 fim ratios are more strongly related to 
the 1.6 fim surface brightness. 



5 DISCUSSION 

The results here have demonstrated that the total stellar popula- 
tion, including evolved stars and not just stars in star forming re- 
gions, plays a significant role in heating dust in nearby galaxies. 
The fraction of observed dust emission that is heated by the to- 
tal stellar population varies as a function of the wavelength. While 
star forming regions are the dominant heating source for dust ob- 
served at < 160 fim, the total stellar population, including evolved 
stars, become the dominant heating source for the dust observed at 
> 250 fim. We have been able to demonstrate that dust is heated by 
the total stellar populations in spite of the fact that the dust emission 
itself tends to be correlated with Ha emission, which should bias 
our results towards finding dust heated by star formation. More- 



Walterbos & Schwerind (1198 



Walterbos & GreenawalM ( Il99i 



over, we have shown that significant dust heating by the total stellar 
population can be seen even in late-type spiral galaxies. 

The results here valid ate the original conclusio ns drawn 
from IRAS data b y iLonsdale Persson & Helou' ('l987j), 

Isauvage & Thuan ( 199^, 

and iKong et a l. (2004) a s wel l 

as t he conclusions draw n from Spitzer data by IHinz et al.i ( l2004h 
and lCalzettietal.1 ( 1201 Oh that evolved stellar populations play an 
increasingly important role in heating dust towards longer infrared 
wavelengths. Additionally, the differences between the results for 
M81, where we were able to demonstrate that the bulge stars were 
heating dust even at 70 fim, and the results for NGC 2403, where 
the dust emitting at 70 fim did not appea r to be influenced at all b y 
the bulge stars, corroborates results from lEngelbracht et al.l ( l2010h . 
who found that bulges influenced dust heating in nearby galaxies. 

Furthermore, the results here corroborate the findings from ra- 
diative transfer models (e.g. Ida Cunha et al.ll2008l : lBianchill20o3; 
Popsecu et al. 2011.) as well as other dust emission models (e.g. 
loraine et al.l 120070 that had suggested that dust absoprtion and 
emission can only be accurately modelled if it contains a dust com- 
ponent heated by star formation and a dust component heated by 
the total stellar population. Both our results and the radiative trans- 
fer analyses suggest that the < 160 fim emission originates from 
dust that is heated by star forming regions located within clouds 
(or clumps) that are optically thick at ultraviolet and blue wave- 
lengths. In these models, only a relatively small fraction of the dust 
is heated by the star forming regions, but the dust appears signifi- 
cantly brighter than the dust outside the centres of the regions, lead- 
ing to higher 70/160 fim surface brightness ratios in these regions. 
Most of the dust mass in the ISM, however, is shielded from the ul- 
traviolet and blue photons from the star forming regions. This dust 
is heated by the diffuse interstellar radiation field, which would in- 
clude light from evolved stars, and has lower temperatures. None 
the less, because a much greater fraction of the dust mass is heated 
this way, the colder dust is the predominant source of emission at 
> 250 fim. Also, if the mean free path of light in the diffuse inter- 
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Figure 16. The observed and reconstructed colour maps for M83. the residual values when the reconstructed maps are subtracted from the observed maps, and 
the Egp I -E-Tota; maps. Each map is 20 X 20 arcmin. See the caption of Figure [Tsl for other information about the layout. 



stellar medium is similar to or smaller than the resolution elements 
in our analysis, the 250/350 /.tm and 350/500 /im ratios can be ex- 
pected to be very well correlated with the stellar surface bright- 
nesses in those resolution elements. The V-band optical depths and 
the dust ver tical scale heights de rived in the analyses of edge-on 
galaxies by IXilouris et"aiT ( ll999h and iBianchil j2007l) imply that 
the mean free paths of V-band photons through the discs of spiral 
galaxies should be roughly between 0.1 and 1 kpc. In comparison, 
the binned data that we used correspond to physical sizes of ~ 500- 



800 pc, and so the 250/350 /im and 350/500 /im ratios should be 
very well correlated with tracers of the total stellar population in 
these binned data, which is exactly what we see. 

While we do produce observational results that are gener- 
ally consistent with many previously-published dust models, we 
do find a inconsistency when comparing our results for M81 and 
NGC 2403 with the dust models for th ese specific galaxies from 
iDraine et al.]j2007l) . lDraine et al.l j2007h did predict from their SED 
titling to M81 that most of the dust emission at > 70 /.tm is from 
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Figure 17. The observed and reconstructed colour maps for NGC 2403, the residual values when the reconstructed maps are subtracted from the observed 
maps, and the Egp / E^otal maps. Each map is 21 X 18 arcmin. See the caption of Figure llSl for other information about the layout. 



a diffuse component, which is largely consistent with the results 
from our surface brightness ratio variations. However, their SED 
fitting to NGC 2403 also predicted that most of the dust emission at 
> 70 nm from that galaxy was from diffuse dust, whereas the sur- 
face brightness variations we observed imply that the star forming 
regions should be the predominant heating source at < 160 fim. 
This demonstrates that, while the general ideas being applied in 
dust models are still valid, further refinement is needed if the mod- 
els are not only going to replicate the overall SEDs of galaxies but 
also replicate colour variations within the galaxies. 

Even though many authors have argued that a significant frac- 
tion of the dust in nearby spiral galaxies is heated by evolved stellar 
populations, other authors have shown that total infrared dust emis- 
sion can be cor related with star formation in nearby spiral galaxies 
jPevereux & Y oung 1990; Dev ereux et al...l995. : .Buat & Xu.l996 : 



iKennicutt et al .120091) . and even some of the early results from Her- 
schel suggest that emission in individual wave bands (i.e. the 100, 
160, and 250 fi m bands) can be correlated with other star formation 
tracers in M33 i lBoauien et alll2010l : IVerlet et aflboid) . a galaxy 
with the same Hubble type as NGC 2403. We think the apparent 
contradictions between the correlations between far-infrared emis- 
sion and star formation found by other authors and the relation be- 
tween dust heating and the total stellar populations that we find can 
be explained in two ways. 

First of all, some of the correlations between far-infrared emis- 
sion and star formation are based on either the far-infrared fluxes 
traced by the IRAS 60-100 /im data or on the total dust emission 
as computed from the Spitzer 24-160 fim bands, both of which 
strongly sample emission from dust on the Wien side of the ther- 
mal dust emission in a regime where we also found that star forma- 
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tion strongly influenced dust heating. Moreover, IRAS-based stud- 
ies simply did not have data at >100 fim wavelengths, which we 
found were dominated by emission from dust heated by the total 
stellar population. 

Second, it is possible that far-infrared emission could be 
linked to star formation indire ctly through the Schmidt law 
( ISchmidl|[l959l ; iKennicuttiri 998bh rather than directly through dust 
heating. Dust emission at far-infrared wavelengths is a function of 
both dust temperature and dust surface density. Since dust should 
implicitly be a tracer of the total surface density of gas in the in- 
terstellar medium, and since the Schmidt law suggests that star for- 
mation is a function of the amount of gas available to fuel it, we 
expect dust emission to be correlated with star formation to some 
degree even if the dust is not directly heated by the star formation. 
The dust colours, however, are strongly affected by dust tempera- 
ture but are independent of dust surface density and therefore will 
depend more on dust heating sources than the emission observed 
in a single wave band. This is probably why we are able to show 
that far-infrared surface brightness in the 250-500 /xm bands is cor- 
related with Ha intensity even though we have demonstrated that 
the 250/350 /im and 350/500 surface brightness ratios are more 
strongly correlated with 1.6 fim surface brightness. Additionally, in 
the case of M81, the 1.6 /im emission is poorly correlated with the 
250, 350, and 500 /xm surface brightness, and yet the 250/350 fim 
and 350/500 ^im ratios are clearly being heated mostly by the stel- 
lar population traced by the 1.6 /im band. In this case the infrared 
surface brightness can only be correlated with the Ha emission if 
the dust surface density is correlated with the star formation sur- 
face density;_Webelieye that this may also be a possible reason 
whv lBoquien et al.l ( l2010t) and lVerlet et al.l ( 1201 O l) can show that the 
100-250 fj,m emission is correlated with Ha and 24 ^m emission, 
although a comparison of the infrared colour variations in M33 to 
tracers of star formation or the total stellar population would be 
needed to confirm this. 

Aside from the use of surface brightness ratios instead of sur- 
face brightnesses or flux densities in our analysis, we think that 
the choice of galaxies that we used in this analysis was critical to 
demonstrating that the total stellar populations (including evolved 
stars) play a significant role in dust heating. M81 and NGC 2403 
are both special in that the stellar emission traces different struc- 
tures than either the dust emission or the star formation. This has 
been critically important for allowing us to disentangle the rela- 
tive contributions from the two different heating sources. M83 is 
a good example of a "typical" spiral galaxy in which the starlight, 
dust emission, and star formation all trace similar structures. Our 
analysis demonstrates that it is difficult to disentangle whether the 
total stellar population or star formation is a greater contributor to 
dust heating in such cases. Other authors working with M83 or sim- 
ilar galaxies will probably encounter similar problems. 

We emphasize t hat we have only foun d that these results ap- 
ply to spiral galaxies. lOalametz et al.l ( [2010h demonstrated that the 
250/500 fim surface brightness ratio depends on 24 fim surface 
brightness in the dwarf irregular galaxy NGC 6822, thus demon- 
strating that the dust emitting at submillimetre wavelengths was 
significantly affected by star forming regions. This is probably a re- 
sult of the galaxy containing relatively few evolved stars that could 
heat dust, although it is conceivable that, in the low metallicity en- 
vironment, the interstellar medium would contain less dust, and so 
diffuse dust would be less shielded from star forming regions. We 
anticipate that similar results would be obtained for other dwarf ir- 
regular galaxies with ongoing star formation. In contrast, we may 
expect that the >70 //m emission seen in elliptical galaxies (when 



the elliptical galaxies produce any far-infrared emission at all) will 
originate from dust predominantly heated by evolved stars, just as 
the dust from the central 3 kpc of M81 appears to be heated pri- 
marily by evolved stars. How our results for spiral galaxies apply 
to luminous infrared galaxies (with total infrared luminosities be- 
tween 10^^ and lO^'^ Lq) and ultraluminous galaxies (with total in- 
frared luminosities above 10^^ L0) is unclear, as the galaxies con- 
tain strong star formation and AGN emission but also have large 
stellar populations and may have large reservoirs of dust that is 
shielded from the star forming and AGN regions. Additional anal- 
yses on the dust heating mechanisms in these classes of galaxies is 
warranted. 



6 CONCLUSIONS 

For M81, M83, and NGC 2403, we have demonstrated that the 
250/350 fim and 350/500 pm surface brightness ratios are more 
strongly correlated with 1.6 /xm emission, which traces the starlight 
from the total stellar populations including evolved stars, than with 
Ha emission from star-forming regions. At shorter wavelengths, 
the total stellar population may continue to be influential in dust 
heating, particularly in early-type spiral galaxy M81, but heating 
by star formation becomes more important, particularly in the late- 
type spiral galaxy NGC 2403. None the less, these results imply 
that, in each of these galaxies as well as other spiral galaxies in 
general, the total stellar populations are a significant if not domi- 
nant heating source for the dust observed at > 160 fim in all spiral 
galaxies, including late-type spiral galaxies. In some cases, virtu- 
ally all of the emission may originate from dust heated by the total 
stellar populations, including evolved stars in the bulges and discs. 

The results here have strong implications for modeling dust 
in spiral galaxies. First of all, dust models need to include sepa- 
rate thermal components for the dust heated by star forming re- 
gions and dust heated by the total stellar populations. When fit- 
ting single modified blackbodies to data, it may be inappropriate 
to force the function to fit 70 fim or shorter wavelength data, par- 
ticularly by varying the index of the power law that describes the 
dust emissivity (see, for example, Dunne et al. (2000), but also note 
the additional analysis in lDunne & Ealesl ( l200ll) ). Models based on 
more complex dust physics need to include components heated by 
red stellar radiation fields that only become predominant at longer 
wavelengths, as is already done in some cases, and these models 
also need to be tuned so that they replicate not only the global 
SEDs of galaxies but also the colour variations within galaxies. 
Template-based SED models need to include an additional com- 
ponent bas ed on dust heated by a quies cent stellar population, as 
is done by Rowan-Robinso n et al. 1(12010'), although these cold dust 
templates should be calibrated using observations of nearby galax- 
ies where the dust emission can be studied on kpc scales rather 
than inferred from observations of unresolved sources. Also, dust 
extinction corrections that rely upon infrared flux to estimate ex- 
tinction correctio ns, as is commonly done for ultraviolet data (e.g. 
iMeurer et alj|l999i) . need to account for dust heating by evolved 
stellar po pulations and not ju st st ar forming regions, a s is also sug- 
gested bv lKong etal] ( l2004h and lCortese et alj ( 1200^ . Finally, we 
are going to suggest caution when using far-infrared emission as 
a star formation tracer for both nearby and more distant galaxies. 
Dust emission in individual wave bands at >160 fim is probably 
related to star formation indirectly through the Schmidt law rather 
than directly through dust heating. Using, for example, 250 ^m flux 
density measurements to determine star formation rates is akin to 
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using CO measurements to determine star formation rates. While 
this may still be an accurate way to measure star formation rates 
under certain circumstances, it will be necessary to understand the 
caveats that affect such analyses. 



ACKNOWLEDGMENTS 

G.J.B. thanks Simone Bianchi and the reviewer for helpful com- 
ments on this paper G.J.B. was funded by the STFC. The reser- 
ach of C.D.W. and K.F. is supported by grants to C.D.W. from 
the Canadian Space Agency and the Natural Sciences and Engi- 
neering Research Council of Canada. PACS has been developed by 
a consortium of institutes led by MPE (Germany) and including 
UVIE (Austria); KU Leuven, CSL, IMEC (Belgium); CEA, LAM 
(France); MPIA (Germany); INAF-IFSI/OAA/OAP/OAT, LENS, 
SISSA (Italy); lAC (Spain). This development has been supported 
by the funding agencies BMVIT (Austria), ESA-PRODEX (Bel- 
gium), CEA/CNES (France), DLR (Germany), ASLTNAF (Italy), 
and CICYT/MCYT (Spain). SPIRE has been developed by a con- 
sortium of institutes led by Cardiff University (UK) and including 
Univ. Lethbridge (Canada); NAOC (China); CEA, LAM (France); 
IFSI, Univ. Padua (Italy); lAC (Spain); Stockholm Observatory 
(Sweden); Imperial College London, RAL, UCL-MSSL, UKATC, 
Univ. Sussex (UK); and Caltech, JPL, NHSC, Univ. Colorado 
(USA). This development has been supported by national fund- 
ing agencies: CSA (Canada); NAOC (China); CEA, CNES, CNRS 
(France); ASI (Italy); MCINN (Spain); SNSB (Sweden); STFC 
(UK); and NASA (USA). HIPE is a joint development by the Her- 
schel Science Ground Segment Consortium, consisting of ESA, the 
NASA Herschel Science Center, and the HIFI, PACS and SPIRE 
consortia. This research has made use of the NASA/IPAC Extra- 
galactic Database (NED) which is operated by the Jet Propulsion 
Laboratory, California Institute of Technology, under contract with 
the National Aeronautics and Space Administration. 



REFERENCES 

Allen R. J., Knapen J. H., Bohlin R., Stecher T. P, 1997, ApJ, 487, 
171 

Arendt R. G. et al., 1998, ApJ, 508, 74 
Bendo G. J. et al., 2007, MNRAS, 380, 1313 
Bendo G. J. et al., 2008, MNRAS, 389, 629 
Bendo G. J. et al., 2010a, MNRAS, 
Bendo G. J. et al., 2010b, A&A, 518, L65 
Bianchi S., 2007, A&A, 471, 765 
Bianchi S., 2008, A&A, 490, 461 

Bohlin R. C, Comett R. H., Hill J. K., Smith A. M., Stecher T. P, 

1983, ApJ, 274, L53 
Boissier S. e al., 2005, ApJ, 619, L83 
Boquien M. et al., 2010, A&A, 518, L70 
Boselli A., Gavazzi G., 2002, A&A, 386, 124 
Buat v., Xu C, 1996, A&A, 306, 61 
Calzetti D. et al., 2005, ApJ, 633, 871 
Calzetti D. et al., 2007, ApJ, 666, 870 
Calzetti D. et al., 2010, ApJ, 714, 1256 
Comte G., 1981, A&AS, 44, 441 

Condon J. J., Condon M. A., Gisler G., Puschell J. J., 1982, ApJ, 
252, 102 

Cortese L., Boselli A., Franzetti P., Decarli R., Gavazzi G., 
Boissier S., Buat V., 2008, MNRAS, 386, 1157 



da Cunha E., Chariot S., Elbaz D., 2008, MNRAS, 388, 1595 

Dale D. A. et al., 2007, ApJ, 655, 863 

Davies J. I. et al., 2010, MNRAS, 409, 102 

de Blok W. J. G., Walter F, Brinks E., Trachternach C, Oh. S.-H., 
Kennicutt R. C, Jr, 2008, AJ, 136, 2648 

de Vaucouleurs G., de Vaucouleurs A., Corwin H. G., Buta R. 
J., Paturel G., Fouque P., 1991, Third Reference Catalogue of 
Bright Galaxies, Springer- Verlag, Berlin 

Devereux N. A., Jacoby G., CiarduUo R., 1995, AJ, 1 10, 1 1 15 

Devereux N. A., Young J. S., 1990, ApJ, 350, L25 

Devereux N. A., Young J. S., 1992, AJ, 103, 1536 

Devereux N. A., Young J. S., 1993, AJ, 106, 948 

Dowell C. D. et al., 2009, in Oschmann J. M. Jr. Clampin M. 
C, MacEwen H. A., eds. Proceedings of the SPIE, Space Tele- 
scopes and Instrumentation 2010: Optical, Infrared, and Mil- 
limeter Wave. SPIE, Bellingham, p. 773136 

Draine B. T. et al., 2007, ApJ, 663, 866 

Drissen L., Roy J.-R., Moffat A. FL,SharaM.M., 1999, AJ, 117, 
1249 

Dunne L., Eales S. A., 2001, MNRAS, 327, 697 

Dunne L., Eales S., Edmunds M., Ivison R., Alexander P., 
Clements D. L., 2000, MNRAS, 315, 115 

Elmegreen D. M., Elmegreen B. G., 1982, 201, 1021 

Elmegreen D. M., Elmegreen B. G., 1987, 314, 3 

Engelbracht C. W., Gordon K. D., Rieke G. H., Werner M. W, 
Dale D. A., Latter W. B., 2005, ApJ, 628, L29 

Engelbracht C. W. et al., 2010, A&A, 518, L56 

Flagey N., Boulanger F, Verstraete L., Miville Deschenes M. A., 
Noriego Cespo A., Reach W. T, 2006, A&A, 453, 969 

Freedman W. L. et al., 2001, ApJ, 553, 47 

Galametz M. et al., 20I0,A&A, 518, L55 

Gamett D. R., Shields G. A., 1987, ApJ, 317, 82 

Gordon K. D. et al., 2005, PASP, 117, 50 

Gordon K. D. et al., 2007, PASP, 119, 1019 

Gordon K. D., Engelbracht C. W., Rieke G. H., Misselt K. A., 
Smith J.-D. T, Kennicutt R. C. Jr, 2008, ApJ, 682, 336 

Griffin M. J. et al., 2009, in Oschmann J. M. Jr. de Graauw M. W. 
M., MacEwen H. A., eds. Proceedings of the SPIE, Space Tele- 
scopes and Instrumentation 2008: Optical, Infrared, and Mil- 
limeter. SPIE, Bellingham, p. 70102 

Griffin M. J. et al., 2010, A&A, 518, L3 

Herschel Space Observatory, 2010a, PACS Observer's Manual, 

Version 2.1, ESA, Noordwijk 
Herschel Space Observatory, 2010b, SPIRE Observer's Manual, 

Version 2.2, ESA, Noordwijk 
Hinz J. L. et al., 2004, ApJS, 154, 259 
Hodge R W., Kennicutt R. C. Jr, 1983, AJ, 88, 296 
Jarrett T. H., Chester T, Cutri R., Schneider S. E., Huchra J. P, 

2003, AJ, 125, 525 
Karachentsev I. D. et al., 2002, A&A, 383, 125 
Kennicutt R. C. Jr, 1998a, ARA&A, 36, 189 
Kennicutt R. C. Jr, 1998b, ApJ, 498, 541 
Kennicutt R. C. Jr et al., 2003, PASR 1 15, 928 
Kennicutt R. C. Jr et al., 2007, ApJ, 671, 333 
Kennicutt R. C. Jr et al., 2009, ApJ, 703, 1672 
Kessler M.F et al. 1996, A&A, 315, L27 

Kong X., Chariot S., Brinchmann J., Fall S. M., 2004, MNRAS, 
349, 769 

Li A., Draine B., 2001, ApJ, 554, 778 
Lonsdale Persson C. J., Helou G., 1987, ApJ, 314, 513 
Lutz D., 2010, PACS photometer point spread function. Version 
l.OLESA, Noordwijk 



© RAS, MNRAS OOQ.fTltlSl 



28 G. J. Bendo et al. 



Markoff S. et al., 2008, ApJ, 681, 905 

Mentuch E. et al., 2009, ApJ, 706, 1020 

Mentuch E., Abraham R. G., Zibetti S., 2010, ApJ, 725, 1971 

Meurer G. R., Heckman T. M., Calzetti D., 1999, ApJ, 521, 64 

Meurer G. R. et al., 2006, ApJS, 165, 307 

Moustakas J., Kennicutt R. C. Jr., Tremonti C. A., Dale D. A., 

Smith J.-D. T., Calzetti D., 2010, ApJS, 190, 233 
Neugebauer G. et al., 1984, ApJ, 278, LI 

Ott S., 2010, in Mizumoto Y., ed., ASP Conf. Ser. Vol. 434, Astro- 
nomical Data Analysis Software and Systems XIX. Astron. Soc. 
Pac, San Francisco, p. 139 

Pilbratt G. et al., 2010, A&A, 518, LI 

Poglitsch A. et al., 2010, A&A, 518, L2 

Popescu C. C., Tuffs R. J., Dopita M. A., Fischera J., Kylafis N. 

D., Madore B. R, 201 1, A&A, 527, A109 
Prescott M. K. M. et al., 2007, ApJ, 668, 182 
Rieke G. H. et al., 2004, ApJS, 154, 25 
Rowan-Robinson M. et al., 2010, MNRAS, 409, 2 
Sauvage M., Thuan T. X., 1992, ApJ, 396, L69 
Schlegel D. J., Finkbeiner D. R, Davis M., 1998, ApJ, 500, 525 
Schmidt M., 1959, ApJ, 129, 243 
Scoville N., Young J. S., 1983, ApJ, 265, 148 
Smith B. J., Hancock M., 2009, AJ, 138, 130 
Stansberry J. A. et al., 2007, PASR 1 19, 1038 
Talbot R. J. Jr., Jensen E. B., Dufour R. J., 1979, ApJ, 229, 91 
Telesco C. M., Harper D. A., 1980, ApJ, 235, 392 
Thim F, Tammann G. A., Dolphin A., Sandage A., Tolstoy E., 

Labhardt L., 2003, ApJ, 590, 256 
Trinchieri G., Fabbiano G., Palumbo G. G. C., 1985, ApJ, 290, 96 
Turner J. L., Ho R T. R, Beck S. C., 1987, Ap, 313, 644 
Verley S. et al., 2010, A&A, 518, L68 
Veron P, Sauvayre A., 1965, AnAp, 28, 698 
Walterbos R. A. M., Greenawalt B., 1996, ApJ, 460, 696 
Walterbos R. A. M., Schwering R B. W., 1987, A&A, 180, 27 
Werner M. W. et al., 2004, ApJS, 154, 1 

Xilouris E. M., Byun Y. I., Kylafis N. D., Paleologou E. V., Papa- 

mastorakis J., 1999, A&A, 344, 868 
Young L. M., Bendo G. J., Lucero D., 2009, AJ, 137, 3053 
Zhu Y.-N., Wu H., Cao C., Li H.-N., 2008, ApJ, 686, 155 



© RAS, MNRAS QOQ.fTHIsl 



